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In the present work, the problem of packetted H.P* 
Turbine Blades is analysed to study its vibrations i.e* to 
calculate the fundamental and a few higher natural frequen- 
cies and to find the corresponding mode shapes. Packetted 
Turbine Blades are made of a number of individual blades 
grouped together by 4 rivetting a steel band at the free 
end. 


Tw modes of vibrations, tangential and coupled 
transverse - torsional modes are considered. The effect 
o.f number of blades in a group and that of band to blade 
stiffness ratio on the natural frequencies of vibrations 
are considered. Finite Element method is used in the 
analysis to obtain consistent mass matrix and corrected 
stiffness matrix. From these Dynamical matrix of the 
system is obtained. Its eigen values and eigen vectors 
give natural frequencies and the mode shapes. The 



stiffness correction matrices are derived to taico intd 
account the effect of centrifugal forces due to the blade 
and the band in both the modes. 

Results are plotted such that the characteristic 
curves give the variation of natural frequencies as a 
function of number of blades .i.n the group for different 
be.nd to blade stiffness ratios. These graphs are to 
enable the designer to avoid a certain frequency of vibra- 
tion b3r proper selection of band properties and / or the 
number of blades to be packotted. 



CHilPTER I 


IFTRODUCTIOK 


It is commonly observed that most of the turbine 
blade failures occur due to fatigue at the blade roots* 
Reliability of operation of steam turbine is of prime impor- 
tance which requires that in design state the mechanical 
reliability of the blade be given a major consideration. 

The fatigue failure of the blade is usually associated with 
the . excessive vibrations in the neighbourhood of its reso- 
nance frequency. For the safe design the designer makes 
sure that the natural f'^^equencies obtained by theoretical 
calculations are reason.ibly away from the forcing frequencies 

To consider tie vibrational analysis of the turbine 
blade 5 some of its important features are presented belowi 

a) Configuration of the turbine blading 

i) Shrouding 
ii) Lacing Wires 

b) Excitations involved in the operation of 
turbinbs 

i) Ihternal 
ii) External 

c) Tunning Process, and 

d) Design Considerations of blade. 
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1.1 Configuration of Blades 

The steam turbines are generally divided func- 
tionally in three parts, namely?- High Pressure (H.P.), 
lilt erinedi ate Pressure (I.P.) and Low Pressure (L.P.) 

The blades are fixed on the rotor at the roots and the 
other end is called the top or free end. The blades used 
in these sections are discussed below. 

1.1.1 High Pressure Section (H.P. Section) 

In the H.P. Section, the specific volume rf 
steam is very low and its expansion is not much. Because 

* 

of these flow conditions, blades used are of short length 
which is of the order of a few inches. The blade is small 
in comparison with the root ~ radius, (1 s 20) giving the 
same tangential velocity throughout its length. From t'^'' 
flow requirements, the entrance and exit angles at all 
points along the blade length are same, and the blade is 
straight and without any twist. 

Shrouding 

For the economical use of the steam, having high 
density, the steam leaking past the blade tip is avoided by 
rivetting a band strip on the free edges of a blade/a group 
of blades, see figure 1.1. This prevention of steam leakage 



3 


is particularly important in the reaction turbines where 
a pressure gradient exists across the moving blades. In 
the high pressure impulse turbines the partial admission is 
employed. The band is beneficial in stiffening the blades. 

The shrouding connecting the blade tops may be 

continuous throughout the periphery of the rotor or in 

segments integral with one or more blades. The shrouding 

is usually attached to the blades by means of tenons cast 

or forged integral with the blades. Tenons fit through the 

holes in the shroud and are rivetted into small spherical 

surfaces. Usually the continuous band is broken into many 

groups of blades with a gap of the order of 0,03'*. This is 

allowed to account for its expansion. The band which runs 

continuously over the periphery, has a tendency to distort 

under the influence of high temperature and may break away 

( 2 ) 

from the tenon causing eventually the blade failure, 

1,1,2 Intermediate Pressure Section Cl.P, Section) 

When the steam passes through I.P. Section the 

steam expands successively. To accommoriate this Increasing 

volume, the steam passage is made wider and wider. This 

essentially requires providing larger and larger blade leng- 

( 2 ) ■ • 

ths along the steam flow direction. The blade height is a 
function of the total annular' area required to pass the 
steam flow, which in turn is a function of mass rate flow, 
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specific volume, and velocity ratio. Due to longer blades 
in the last ring of the I.P. Section the blades are slightly 
twisted. The blades still need shrouding.. The lengths of 
the blades are of the order of 4 to 6 inches. 

1.1.3 Low Pressure Section (L.P. Section) 

In this portion the specific volume becomes very 
large. The blade lengths are of the order of 10 to 12 in- 
ches. The blades become appreciably long and as such entra- 
nce and exit blade angles vary along the blade lengths. 

For this it requires that the blade should be twisted. The 
blades are also made tapered by cutting down the blade mass 
at the ends so that the level of stresses due to centrifugal 
forces at the root of the blades is reduced. 

Lashing Wires 

One of the reasons for using shroud is to hold 
blades in the correct position even under the influence of 
steam forces. This is also achieved by Lacing or Lashing 
Wires. Shrouding introduces high centrifugal forces where 
as lacing achieves the purpose of alignment without appre- 
ciable addition to centrifugal forces. 

Formerly lashing wires were strung through the 
bossed holes, machined in the blades. The present practice 
is to weld short pieces of lashing wires to each side of the 
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blade. When the blades are machine - finished and stresses 
are relieved, these lashing lATires touch each other and are 
joined either by a weld or by a welded sleeve. The main 
disadvantage of lashing wires is that it disturbs the flow 
pattern and contributes to the vibration tendencies. To 
avoid these, ao - rofi - ll section is used for the lashing wires. 
Lashing wires could be provided at two or three places along 
the length of the blade in order to give it sufficient 
rigidity. 

1,2 Excitation Forces Developed in the Turbines 

Having presented the configuration of the blades, 
we present below different types of disturbing forces deve- 
loped in the steam turbine. The fundamental forcing frequen- 
cy imposed on the blade is equal to the turbine rotor speed 
where as its higher harmonic excitation frequencies are 
due to varying steam flow from nozzle to nozsle. However 
these exciting frequencies are not necessarily same for 
different stages of the turbine. Different kinds of exci- 
ting forces are given below. 

1.2,1 Internal Excitations! 

In stationary steam turbines the main source of 
excitation is the steam itself. Even in complete circumfe- 
rencial steam admission, the steam flow is not uniform due 
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( 3 ) 

to the following practical structural features 


(a) 


The runner blades experience steam flow of vary- 
ing velocity as these pass across a stationary 
guide vane , having maximum velocity in the central 
zone of the stream and minimum velocity at the 
portion issuing out from the neighbourhood of the 
walls of the nozzle, see figure 1.2, This velo- 
city variation imposes a forcing frequency given 
l3y, 



c.p.s 


( 1 . 1 ) 


where n = no. of nozzles for a ring 

N = speed of the rotor in r.p.m. 
f = forcing frequency 

(b) Several discrete steam inlet pipes 

(c) Discrete steam extraction pipes 

(d) Lack of uniformity in the nozzles or guide blades 
may also become a source of uneven forces. Due to 
discontinuity at the horizontal ooint, it is diffi- 
cult to keep the nozzles at the same spacing. 
Because of the changed positions of nozzles at the 
horizontal joint, the runner blades experience addi 
tional disturbing forces twice per revolution. 


(e) Structural members in the inlet and exhaust pipes 
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(f) Manufacturing variation in the nozzles and blades 

(g) Impulse Excitation 

this occurs in the H.P, stage due to partial 

admission of steam and the blades experience 

(4) 

sudden application of the full force, 

1,2*2 External Excitations 

There are also some forces external to the rotor 
which cause excitation and these are as. follows : 

(a) Due to lack of cyliMrical nature of the journal* 
This causes a motion which may set blades in vi- 
brations. 

(b) In marine applications, the turbine - gear - 
propeller system develops strong torsional osci- 
llations of low frequency which are imposed on the 
blades. 

(c) Due to gear tooth error, exciting forces of higher 
frequencies are also set up. 

The excitations mentioned in (b) and (c) are not 
present in the stationary turbines due to direct coupling of 
turbine rotor to the generator. 
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1.3 Tunning Process 

Tunning process is a process of varying physical 
parameters of a "blade like length, breadth, thickness and 
its taper. This may be done by using a parameter or a com- 
bination of parameters subject to certain constraints. 
Tunning is necessary in order to set a particular natural 
frequency. This is carried out after the blade has been 
designed from the flow and strength considerations. In the 
case of banded blades this can be achieved by additional 
means, such as varying the geometrical properties of the 
shrouding or varying the number of blades being grouped. 

The tunning process is conveniently carried out 

by- means of Campbell diagram named after late Wilfred 

(5) 

Campbell of G,E,C, 

In the Campbell diagram^^’*^’^^ important harmonic 
frequencies, (usually cycles per minute) are plotted against 
turbine revolutions per minute, see figure. 1,3, These har- 
monics are represented by straight lines passing through 
the origin with their slope equal to their ith order of 
harmonic. In addition to these the natural frequencies of 
the blades fj_ are plotted on the same graph. The inter- 
sections of these two sets of lines fj_ and h^ indicate the 


resonance* 
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With the help of this diagram, the designer avoids 
the coincidence of important natural frequencies f^, (funda- 
mental or higher) with the forcing frequencies of excitatiofis. 

/ 

Ideally, the intersecting point a^ of the hamonic h^ and •" 
the rotor- speed N should lie symmetrically ..about the inter^ 

V Sections b^ and b^^^. of natural frequency ,f^ and the rotor 
speed ordinate M. . 

' ( I 

Vibration per revolution, , denoted by VR, is a 

term defined as the ratio of blade frequency f to the rota^ 

■ i 

tional speed of the rotor N, This may or may not |)e an •, 

t \ 

integer. Vibrations per revolution VR of 3 should be avoi-* 
ded because it means that the forced frequency has its third 
harmonic hg coinciding with one of the fj _5 similarly VR = 4.9 
means that fifth harmonic h 5 is very near to . one of the na- 
tural frequencies of the blade and should also be avoided. 
Vibrations per revolution of 2,5 means that the intersection 
points ag and ag are equally away from b^ and b^^^. Hence 
the values of VR 2,5, 7.5 etc. should be ideally achieved* 

For short blades centrifugal action has negligible 
effect on natural frequencies f^^s. Whereas for long blades 
it is prominent. Thus the natural frequency of the blade is 
independent of N for short blades and depends upon N for 
long blades. In the tunning process the straight lines/ 
curves for are either raised or lowered so as to obtain 
the ideal VR. 
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1,4 Design Considerations of Blade 

The blade is primarily designed from the flow and 
strength considerations. The forces acting on the blade are 
of two types namely static and dynamic loading, 

1.4.1 Static Loading 

The static loading on the blade is due to steam 
pressure and the centrifugal force. In general, these forces 
are high but fortunately these are of static nature and do 
not cause difficulty in designing for the same. 

1.4.2 Fatigue Loading 

Several turbine accidents have been directly attri- 
buted to blade vibrations caused by synchronism of certain 
periodic disturbing forces with the natural frequency of the 
blade. Hence one must be particular to take into account 
these vibratory forces to avoid the possibility of fatigue 
failure. 

As mentioned before, the blades are subjected to 
the vibrations caused by irregularities in the fluid flow, 
non- symmetry of stationary guide vanes, radial pressure dis- 
tribution, shock and partial admission, etc, A similar 
effect is also produced in the gas turbines due to seperate 
combustion chambers. 



Lashing wires or shrouds add stiffness to the bla- 
des and in some measure these can be used to change the 
natural frequency of the group. ¥e have seen that the blade 
configuration varies from stage to stage in a section. Thus 
blades in different rings have different sets of fundamental 
and higher natural frequencies. Similarly the number of 
nozzles per ring may vary from stage to stage and here too 
we get a spectrum of forcing frequencies. It follows that 
for a given rotor speed, blades jn. one or more than one ring 
may run with natural frequencies equal to or in the neigh- 
bourhood of the forcing frequency. In case of constant 
speed stationary turbines, it is possible to tune blades of 
all the stages. 

It is difficult to achieve the same in case of 
variable speed marine turbines. For these, tunning process 
is carried out for the maximum operating rotor speed so as 
to keep the vibratory forces at the minimum even though the 
centrifugal forces are high. For any other operating speed, 
the resultant forces are within the designed limit. 



CHAPTER II 


REVIEW OF THE PREVIOUS WORK 

Turbines are generally divided functionally in /' 
three different parts namely H.P. , I.P. and L.P. Sections. 

In the previous chapter we have mentioned about the types 
of blades used at different stages of the turbine. Blades 
used for H.P. section are short and straight and are always 
bonded. Blades in I.P. section are of medium heights, some- 
times slightly twisted, and are usually banded. The blades 
in the L.P. section are long, tapered and twisted. Generally 
these blades are held together by lashing wires. We present 
here a short review of the work done in the area of vibra- 
tion of turbine blades, with particular emphasis on tapered 
twisted blades and banded H.P, blades, 

( 6 ) 

William Carnegie and J, Thompson considered the 
effect of width taper on frequencies and mode shapes of 
straight rectangular blade . This cross - section is selec- 
ted to isolate the effect of taper from those of parameteitJ 
like twist and nonsymmetry of section. Width taper is 
varied ov$r a wide range , from - 0,07 to + o.l inch per 
inch of blade length. Finite difference method is used for 
the derivatives which reduce the equations of motion to a 
set of algebraic simultaneous linear equations which forms 
an. eigen value problem. The results obtained from the eigen 
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value problem are compa.rable with the experimental results. 
The authors concluded that the natural frequencies decrease 
with the increasing positive width for all inodes of vibra- 
tion and the change in frequencies is more sensitive to 
negative range of values. 

(7) 

Rosard D.D considered the natural frequency of 
twisted blade of rectangular cross - section v/ith particular 
interest to investigate the effect of width to thiclmess 
ratio. In general, when two systems of nearly equal natural 
frequencies 00 ^ '><■02 are coupled, the natural 

frequencies ^c2 ^^cl^ coupled system are 

such that ^g 2')'‘^2 * case of beam 

vibration problem the two vibrating systems involved are the 
second beam mode vibration in the flexible direction and the 
first mode vibration in the stiff direction. The coupling 
is introduced by the twist of the beam. The author found 
that the coupling becomes appreciable for the width to thick- 
ness ratio in the range of 4 to 12, The problem is analy- 
tically solved by the lumped mass system and by applying 

( 8 ) 

Myklestad Method , The values obtained are compared with 
the experimental results and the errors involved are in the 
range of 0.1 to 2 per cent. 

(Q) 

B. Dawson has used Rayleigh - Ritz method to 
calculate the natural frequencies of coupled bending - 
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bending vibrations of pre-twisted cantilever beam* The beam 
of aerofoil section is considered. The expressions for 
kinetic and potential energies are derived and are as 
follows, 

( 'Si H El "1 

PE = V = j ^ (y")^ +E y" z" + — ^ (z”)^| <ax (2.1) 

L 

KE = T = ( (y^ + z^) dx (2,2) 

The deflection curves y (x), z (x) are assumed in a way to 
satisfy the geometric boundary conditions and need not satis- 
fy the static boundary conditions. The boundary conditions 
used are as follows s 

1. Geometric B,C . 

at X = O 5 z = 0, y = 0, z' =0, y' = 0 , 

2. Static B.C. 

at X = L 5 z” = 0, y" = 0, z'” = 0, y” ’ = 0 

(2.3) 

The torsional frequency of straight unifomn beam, 
encastre at one end and free at the other end, allowing for 
Saint Venant resistance to shear only, is given by, 
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whe re n 

Cl 

L 


odd Integer 

stiffness constant 
3 

G a b / 3, for rectangular section 
length the bar 

polar mass moment of inertia 


Equation (2.4) is modified by W. Carnegie^^^^ using 
an energy method, and has produced a correction factor which 
allows for the bending of longitudinal fibres during torsion. 
The frequency equation is given by 


f 


n 





^ 

Cl 


(2.5) 


The corresponding equation for pre-twisted beams 
of uniform cross section is obtained by replacing^ Cl by 

Cl + G3 (d^/ dx)^ (2.6) 

where pi - angle of pre-twist ■ 

= (x/L) 

being the total pre-twist. 

W. Carnegie^^^^ developed total energy expression 
for a thin blade vibrating in coupled bending - bending - 
torsion giving fibres due to torsion* Use is made of ex- 
pressions for PE and KE, equations (2.1) and (2,2), for 
bending - bending of the blade, derived in (9) and (10). 
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For calculating expressions for total potential anu kinetic 
energy and for bending - bending torsion of the thin 
blade , additional terms are 


V, 


' _Q_1£L1 - E (en(Q»') /^2 3 

2 m 


b t db 
A 


dx 


^ Ic.G 


o. 


w 


dx 


(2.7) 


( 2 . 8 ) 


For aerofoil' sections, the centre of twist is away from its 

C.G, and let the location of shear centre be at distances r. 

r in the y and z directions respectively. Then the expre- 
z 

ssion for total K.E, becomes 


y’ 


T = 


( X + 0) 

r y 


2g 


(y + r Q)2 1^^^ (Q)‘ 

. + . * 


X 


2g 


2g 


dx 


(2.9) 

Expression for Lagragian function L is obtained. The . 
Hamilton’s Principle is used for 'getting equations of motion 
and the boundary conditions. The equations of motion are 
given by 


dx 


E l - E 1 v" 

a Xyy 4 XK ^ 


w z ^ TiT ^y p ^Q 
g ” S 


(2.10) 


r B y" - E 1,, z" 


dx 
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A special case of straight blade, of uniform cross 
section having symmetry about both the principal axes, is 
considered. The torsional vibration frequencies are cal- 
culated which agree closely with those obtained by Rayleigh 
method and experimental method. 


( 12 ) 


A.P, Duggan and H.A. Slype r 
finite difference method to equations of motion. 


used modified 

( 11 ) 


ro ^4. 2d C2 d^© 

, 4 dx"' .^3 

dx dx 


+ f dC2 . 03 . 01 

L dx^ ^ 


d © _ d Cl d© 

Jdx^ dx * dx 


d G3 (doC ^= i60^q 
dx dx dx p 


(2.12a) 


and computed the natural torsional frequencies of vibration 
of the blade. The blade of uniform cross section is con- 
sidered, Allowance is made in the analysis for the bending 
effect of longitudinal fibres during torsion. The results 
are compared with those obtained experimentally which are 
within accuracy of 1^. The error is minimum when the beam 
is divided into elements varying from 30 to 40, The results 
obtained for rectangular section are comparable with those 
obtained by Carnegie and the exact solution. 

" ('1 3'\ ■' 

W. Carnegie derived general equation of motion 
for bending - bending - torsional vibration of pre-twisted 
cantilever beam. The analysis takes into account the 
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effect of rotary inertia, shear determination and the effect 
of tending longitudinal fibre due to torsion. The expre- 
ssion for P.S and K.E are derived. The variation principle 
is used to derive equations of motion and the boundary con- 
ditions. Two important assumptions made 

1, Strain energies due to bending, torsion and shear 
are independent of each other. And, 

2, Plane cross-sections remain plane before and after 
deflection and no warping is produced. 

(14) 

J.S. Rao and W. Carnegie applied Galerkins 

procedure for finding the natural frequency of straight 

imi form blade with assymmetric cross-section executing 

coupled bending - bending - torsional vibrations. The 

results obtained are comparable with those obtained by a 

(15) 

numerical procedure developed by Rao and also with 

( 16) 

the experimental results obtained by Carnegie & Dawson . 

(17) 

¥. Carnegie also developed equations of motion 
of a rotating blade executing small vibrations.' In the 
paper, expression for work done against centrifugal forces 
is derived. Expressions for total P.S. and K.E. are ob- 
tained. Variation principle is used to get characteristic 
equation of motion for rotating system of general aetofoil 
section. Rayleigh’s method is used to find the fundamental 
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frequency for the case of uniform hlade of symmetrical 
cross - section. 

( 18 ) 

W. Carnegie in another paper derived diffe- 

(17) 

rential equations of motion of the same system by in- 
cluding the effect of rotary inertia and the correction 
factor due to bending of longitudinal fibre during torsion. 
The author found that the natural frequency depends upon 
stagger angle 0 and that the fundamental frequency is 
unaffected by 

1) pre -twist and 

2) coupling between bending and torsion. 

(16) 

W, Carnegie and B, Dawson considered the 
vibration of straight blade of aerofoil cross-section. The 
theoretical procedure essentially consists of transforming 
equations (2.10), (2,11) and (2.12) to a set of 10 first 
order differential equations. Solution of this set is 
obtained by Runge Kutta step by step numerical integration 
procedure. For the case of, straight beam equations (2.10), 
(2,11), (2,12) are reduced to 


E I 


XX 


d^y 

d? 


g 


2 

w 2 n 9 


(2.13) 
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7«r^y p -z ^ ^z P^ y 


(2,15) 
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\ X 

Solution of tliG form z = A e is assumed. On 

substitution In. the set of ten ordinary differential equa- 

2 

tions, a soco.ndary quintic equation is obtained in A . 
Using ten boundary conditions, the set of simultaneous 
equations are solved for constants of integration. 


The results obtained by Runge Kutta method are 

comparable with the analytical and the experimental results, 

obtained by Carnegie*'^^^ and are within 0.25?^ accuracy. 

( 2 ) 

G, Esakson and J.G. Eisley ^ discussed both the rotating 
and non rotating twisted blades for finding the natural 
frequencies .in bending. Both cantilevered and articulated 
blades are considered. Offset of the support from the axis 
of rotation has also been included, A relation developed by 
Lo and Renbarger^^^^ for the effect of blade root angle 0 
on the natural frequencies is found to provide useful appro- 
ximation. The analysis makes use of method developed by 
( 21 ') 

Targoff which is an extension and adoptation of Holzer- 
Myklestad method. The Rayleigh Southwell approximation, 
according to, which the PE of the system is found by intro- 
ducing the work done against the centrifugal force, is used. 
The relation between the frequencies of rotating and non- 
rotating blades is given in the form. 






( 2 . 16 ) 
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whe re 


^rn 

PNn 


k 


= natural frequency for rotating blade, 

*1* To 

= n natural frequency for non- rotating blade. 
= A constant for n"'^^ mode and is given by 

!f rX 




i 2 ' 

I f X djz ( yj^) + ( z^) 


- oy 




n 


dx 


n 


1 j 


I ^ 2 ^ 2 . , 

I t' ( ^ ^ ^ 


(2.17) 


"n 


A case study of overhead helicopter propeller blade is 
presented. 


The author concludes that, 

1, Fundamental natural frequency of cantilever and non - 
rotating articulated blade is almost i nde pendent of 

. twist if the ratio of raaDor to minor stiffness ratio 
exceeds three. 

2, The second and the third natural frequencies of non - 
rotating cantilever blade decreases with increasing 
twist if its maximum stiffness is very large compared 
to its minimum’ stiffness. 

3, Fundamental natural frequency of a rotating cantilever 
beam depends primarily on the root blade angle and 
secondarily on the twist and 

4, Offset of the root support increases the frequency. 
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The theoretical analysis of the tvristed turbine 

blade 5 both rotating and non rotating, was undertaken by 

( 22 ) 

R.C . Diprima and G.H. Handleman . The co-ordinate sys- 
tem used consists of a triad of orthogonal unit rectors, 
fixed at the root of the blade and another moving triad 
sliding along the blade length vj'ith angular speed equal to 
the twist per unit length* Effects due to rotary inertia 
and the shear deformation are excluded. The cantilever of 
a general cross-section is considered. The beam is assui^od 
to undergo small transverse vibrations. Quantities like 
deflection, strain and stress are expressed as vector quan- 
tities. Equations of motion are obtained by equating the 
restoring forces to the inertia forces. The same equations 
of motion are also derived by using variational principle. 
The solution is obtained by assuming a set of fourth order 
polynomials for the two deflections and finding the cons- 
tants of polynomials by trial and error method. The case 
of a propeller blade is considered for large width to thick- 
ness ratio-, The results obtained agree closely with the 

( 23 ) 

experimental results . 

So far we have reviewed the work done for a 
single rotating or non- rotating twisted cantilever blade. 

Now an account of work done on the packetted turbine blades 
is given. Only a few papers have appeared in the litera- 
ture and these concern mainly non-rotating banded blades. 
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( 24 -) 

M*A. Prchl has attempted the problem of 

packetted turbine buckets for the H.P, Section wherein the 

bunded group consists of evenly spaced identical, parallel 

buckets. For generality blade with hinge flexibility at 

the root and non uniform cross section is considered. 

Longitudinal displacements of the blade are assumed negli- 

gible with respect to their transverse deflections. Shear 

deformations and the rotary inertia effects are discarded, 

( 26 ) 

His analysis follows the approach used by Smith . The 
principal axes of moment of inertia of the cross-section are 
assumed to be parallel and perpendicular to the plane of 
the wheel on which the banded group is mounted. The centre 
of txfist and its centre of gravity are assumed to coincide, 

A modified Holzer technique is applied for the two modes 
of vibration encountered in the bonded blades. These are 
as follows , 

1. Pure tangential mode wherein all displacements are 
parallel to the plane of rotation and 

2, Coupled transverse and torsional mode in which the 
displacements are perpendicular to the plane of disc 
and the angular motion is about the longitudinal axis 
of the blade. The coupling between the transverse 
and the rotational displacement is due to the band 
connecting. the adjacent blade tips. The use is made 
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of symmetry and anti symmetry of the banded blades for 

groups having odd and even number of blades. This helps 

in reducing the order of frequency determinant to almost 

(25) 

half. In his next paper , he presented a case study. 
The paper includes a procedure to determine vibration 
amplitude and stresses at resonance. Damping is assumed 
to be small. To determine the 'stress level, the energy 
input to the bonded group from a prescribed form of nozzle 
stimulus is equated to the energy dissipation in damping. 


(5 ) 

Wilfred Campbell did the experimental inves- 
tigations for L.P. blades grouped by lashing wires and found 
the effect of speed of rotation on the natural frequency. 

An Imperial relation between p^^^ and p^^^ was presented 


2 _ 2 , -n 

p = P.T + B K 
^rn -^Nn 

where B is obtained from the test data for the bonded blades. 
B for different cases is presented below 


(1) B = § a ^ § (2.18) 

for a bar of miform section and 

(2) B = + I (2.19) 

for a rectangular bar of constant width and 
tapering to zero in thickness at the free end. 
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The paper also considers the tiinning process in detail* 
W.C, Heckman in part II of the same paper presented the 
effect of number of buckets tied together by lacing wires 
for L,P. section. The plotted results show that the natu- 
ral frequency increases duo to banding. 



CHAPTER 3 


THEORETICAL ANALYSIS 

In the previous chapters, we have stated that the 
blades in the H.P, section are short, straight, irntwisted 
and invariably banded. The blades have aerofoil cross 
section and the shear centre of the cross - section is 
slightly away from the C.G. The principal axes of moment 
of inertia are not exactly parallel and perpendicular to 
the plane at the section. In actual practice the band and 
the blade tip joint may not be perfectly rigid. 

As mentioned in the review, the problem of banded 

(24) 

H.P, blades was solved by M.A. Prohl by Holzer' s Method, 
He did not consider the effect of centrifugal forces. In 
the present work, the problem of banded H,P. blades is being 
analysed for vibrational analysis by using Finite Element 
Method, incorporating in the effect of centrifugal forces. 
The object is to collect necessary data for design purposes 
of banded blades. For this purpose natural frequencies and 
mode shapes of banded blades are obtained for 

a) number of blades in the group being banded 

b) band to blade stiffness ratio 

c) the rotor speed, • 
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3,1 Assiomptions 

In the H.P, section, the blades of different 
stages are fixed to a common rotor and do not transmit 
vibrations to the rotor^^^. This implies that the vibra- 
tions of the blade/bnnded blades is an independent phenomenon. 

The H.P. blades described above, are idealized 
for carrying out analysis and the assumptions made are 
mentioned below. 

1. The blades banded together have the same configu- 
ration and physical dimensions. 

2, For the blades, offset of shear centre from the 
C.G. is very small and are assumed to coincide* 

The cross section is assumed as rectangular having 
the same inertial properties as those of original. 

3, The staggering angle of the blade root being very 
small, the principal axes are assumed to be para- 
llel and perpendicular to the plane of wheel or 
the plane of rotation, 

4. The blades are assumed to be rigidly fixed at the 
root and no hinge flexibility is considered; as 
such the blades are treated as cantilevers. 
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5* The ri vetted ;ioint between the band and the blade 
tips is assumed to be perfect. This does not per- 
mit relative motion between the band and the blade 
at the joint. 

6. The radius of rotor at the blade root is very large 
compared to the length of H.P, blade (nearly of the 
order of 20 : 1). This implies that the banded 
blades; in a group of 6 to 8, can be considered 

to be parallel to each other. This facilitates 
in choosing the nodal co-ordinate axes parallel 
to each other. Thus no transf ormat3.on is needed 
for the superposition of individual elemental 
matrices. 

7. Shear deformation and rotary inertia effects are 
not considered. 

8, The blades being of short length, the longitudinal 
deflection is assumed to be negligible in compa- 
rison to transverse and tangential displacements. 
This helps in reducing the order of stiffness and 
inertia matrices for the whole system and 

9, External forces considered in the analysis are 
those which arise from centrifugal action and 
are of two types. First one is the body force 



32 


distributed along the blade length. The second one 
is force which is produced due to centrifugal force 
of the mass of the band, concentrated at the blade 
tip. 

3.2 Problem Formulation 

Finite Element Method is used to formulate the 
equations of motion. In the banded blades, the blades and 
the band are divided Into finite elements. A certain co- 
ordinate system is defined at each nodal point depending 
upon the mode of vibration under consideration. For each 
such element, the equation of motion^^"^^ is given by 



a U dv + 


b*^ U dv 


v-" 


P + 


0 ds + 


X dv 


V 


(3.1) 


where 

- mass per unit length 

a - a square matrix which relate s the local co-- 
ordinates (x and Uy) in terms of nodal dis- 
U (U2, Ue, US, U12) 

b - a square matrix relating displacement and 
strain vectors 

P - External concentrated force 
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xb - a square matrix relating stresses and 
strains. 

0 »- surface forces, which are function of local 

co-ordinates 

X - a vector of external body forces expressed 
in terms of local co-ordinates. 

U - vector of nodal displacements. 

U - second derivative of nodal displacement U 
with respect to time. 


The equation (3,1) for i"*^^ element of the beam becomes 


«1 ^1 Ki = P3_ + 




^ 0^ ds + 


faT 




1 X, dv 
1 


where M. = 
1 


K. = 
1 




V. 


a . dv 
1 1 


b. (x ) . dv 
i b 1 


> by definition 


(3.2) 

(3.3) 

(3.4) 


and are called inertia and stiffness matrices. Similarly 


(P ). 
s 1 


Si, 


T 
i 


a. 0, ds, (P, ) 


b i 


X, dv 


are called as generalized surface and body forces respec- 
tively. The integrations are carried o ut for the beam 
element (Fig”. 3.1-b) having six degrees of freedom at each 
node. These two matrices for the beam element are of the 
order of 12 x 12 and are glv^en in the Table 1 & 2 respectively, 
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STIPFHESS MATE3X FOR M EIiEMERT OP BEM WITH 6 DEGREES CP PEEEDCM 

(Ref. Pig..3,1-b) 


M 

1 


0 


12EI 


0 

0 

0 


12EI 

0 


0 


0 


-6EI 

0 


GJ 

1 

0 


"M 

1 


6EI 
z 

0 

-12EI 


0 0 


4EI 


0 


0 0 


0 


-12EI 6EI 

0 


4EI 

z 

1 

0 

-6EI 


1 


0 


12EI 


r 


SIMM, 


12EI 


0 


0 


0 


0 


6EI 


•^GJ 

1 


-6EI 

— ^ 0 
1 


2EI 


0 


2EI 


0 0 0 


0 0 

0\ 0 
-6EI 


0 


0 


6EI 


1 


0 


4BI 

0 


0 0 


4BI. 



TJfflIE - 2 


MASS MATHH FOR Ml EIEMEIT OP BEAM TiTITH 6 lEGEEES OP PEEBDCM 

(Ref. Pig, -3.1-c) 



HO 
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M is called consistent mass matrjoc. Unlike the 
mass matrjjc in lumped mass system, M obtained from equation 
(3,3) is not a diagonal matriz. Having obtained IL , for 
all i‘s superimpose these in an appropriate manner to ob- 
tain mass and stiffness matrices M and K respectively for 
the whole system. Then the equations of motion take the 
f o rm 


M U + KU = P + P +P, (3,6) 

s b 

where U defines all the nodal displacements- and M, K, P, 

Pg and P-j^ correspond to the coordinate system of U, Ijfcen 
only the body forces due to centrifugal action are consi- 
dered, the equation (3,5) becomes 

M U + K U = P t P^ 

= P + U (3.6) 

or M U -4- (K - K^) U = P 

M ii + K U = P 

where K is the corrected stiffness matrix and K = K - 
is the stiffness correction matrix. For P =0, let the 
harmonic solution be of natural frequency p, then equation 
(3.6) reduces to, (K is represented as K without any loss 
of generality) 

-p^MU-<-KU = 0 

V ^ ^1: ^ : '3 

Premultiplying by K throughout and rearranging, wo get 
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or 

DU = X U (3.7) 

where 

D ' = K~^ M : 

\ ~ —1—. (3.8) 

P 

Thus the prohlem of the vibration of the system reduces to 
■an eigen value problem, equation (3,7). Matrix D whose 
eigen values are to be found is called Dynamical Matrix and | 

is given by the equation (3,7), The eigen values give 
the characteristic values of the system and these give its 
natural frequencies. The eigen vector U corresponding to mg 

gi-ves the mode shape of vibration for the natural frequency 
given by \ . I 


Fig. 3..2-a shows the configuration of the turbine 


blade mounted on the turbine rotor. The orthogonal co-ordi- 
nate system (X, Y, Z) are selected as sho\m in the figure. 
T%\ro modes of vibrations are considered. Fig^ 3.2-b shows 
the co-ordinate system chosen for the tangential mode of 


vibrations wherein the displacements (Ug, UgjUgjU^g) occur 
in the X - Y plane, Figi,. 3.2-c shows the co-ordinate system 


for the coupled transverse and the torsional vibrations 
wherein the displacements (U 3 , U^, Ug, Ug, U^^) occur, 

in the X - Z plane and the angular rotations about the indi- 


vidual blade axis. 
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3.3 Derivation of Stiffness Correction Matrix 

In the last section, we have presented how the 
mass and stiffness matrices for the system are found from 
equations (3,3) and (3,4). It is also mentioned that depend- 
ing upon the external forces, equation (3,5) gets modified. 

In this section we will consider tho effect of external 
forces in two modes of vibrations namely tangential and 
coupled transverse - torsional modes. Equation (3,5) is 
reproduced hero as 

M U + K U = P + X dv 

= P + (3,9) 

Effects due to concentrated force P and the distri- 
buted force X are considered soperatoly and these are found 
in both tho modes of vibration. Thus the four correction 
matrices derived are as follows, 

a. Tangential mode 

- takes into account , the effect of distributed 
body force of the blade element due to centri- 

Kg - takes into account the effect of rotating mass 
of the band connected at the end of the blade. 
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b. Coupled transverse and torsional mode 

Kg - takes into account the distributed body force of 
the blade element due to rotation, 

K4 - considers the effect of rotating mass of the 
band connected at the end of the blade. 

These are derived below. 


3,3.1 Expression for K 




1 
3-3q 


Fig', shoxiTs an element of the blade marinted 

on the rotor and executing tangential mode of vibration in 
the X-Y plane. The rotor rotates at a constant angular 

velocity of jCl , The nodal co-ordinates U„, U^, U . U 

y2\, 

have already been shown in figure 4^). The inner end of 
the oloment is at a dj.stance R from the axis of rotor* 


For the deflection curve corresponding to tangen- 
tial 'vibrations ‘the distribution of body force X is found 
boloxf t 
Let 

dF = Force in the radial direction OC duo to a 
small element dx • situated at a distance x 
from the inner end A. 

dF^ = Component of the centrifugal force dF along 
X axis 
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dF = com-oonont of tho centrifugal force dF along 

y 

Y axis. 

U = Tangential displacement in the Y direction of 
y 

the oleEient situated at distance x. 
m = mass per unit length 

A = Cross-sectional area of the Llado (which is 
constant in this case). 

When the rotating hlade is not vibrating, the 
body force acts along the blade axis x but when the blade 
starts vibrating, the transverse displacement gives rise 
to component dF^ in the Y direction. This vertical compo- 
nent of centrifugal force affects the vibrational phenomena 
by introducing additional restoring forces which is a linear 
function of displacement U .. For an element at a distance r 

V 

and having a transverse displacement U^, 

/ 2 2 ^ 

the radial distance OC r + 

Hence tho centrifugal force dF, its components dF^ and dF^ 
and the body couple are 

dF = 

dF^j - m j\. 

dM = 0 

xy 


X 




= m r -11 




m U 2 

y 
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The external hody force 

-L * 


d F. 


X 


d F 


IS 


d M 


'xy 


M r-A. 


X = <1 M U 

1 1 y 

0 


2 > 


(3*10) 


U Is expressed in nodal displacements (Uo , U , U ) 
y ^ 2’ 6’ 8’ 12 


where , 

2 3 — 

U = (1 - 3^ + Sf ) U + cf - 2 1 +f'')lU_ 

y ' 2 ^ o 


2 ^.3 

y . <2' 


4- (3 I 2 Ug + (- 1 

where 1 = length of the finite element AB 

f = x/1, X being distance measured from A. 


(3.11) 


When the expression (3,10) for X^ is used in equation (3*9), 
it becomes, 


M U K U = 


a X dv 


V' 


1 ’ 


say 


(3.13) 


where the integration is carried over the volume of the 
blade. When the value of U from expression (3.11) is subs- 

y ■ ■' : 

tituted in the E.H.S. of eauatioxi (3,12) , F^* Is found to 
consists of two parts, namely . 

F^ = component of F^* which is a function of nodal 
displacements, Ug, Ug, Ug, 

F' = a constant vector 


F^* is then written as, 
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F * = K IT + F’ 

1 1 1 

where = squa,rG matrijc. 

Replacing F^* in equation (3.12) 


(3.12 ) 5 we ge t 


by the R.H.S. 


(3,12) 


of 


M U + K U = U + F'^ 
or 

# 

M U + K. U = F^ ’ 

IL 

whe re 

= I^K - K^jand is called the corrected stiffness 
matrix 

K = original stiffness matrix for stationary element 
= stiffness correction matrix. 

Tho correction stiffness matrix is given in table 3a. 


3,3.2 Expression for Kg 

Fig". 3,3-b' shows an ole me nt at the tip of a blade. 
The force F, due to centrifugal effect of rotating mass M 
of the band, is shovra to act at tho tip B. ¥e find below 
the expression for Kg due to the body force of tho band 
which acts concentrated on the blade end. Hence and P-j^ 


are zero 
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a) STIPPNESS CORHECTIOW MATRIX 

* 1 

I A b| 


Where , 


1^ 

(4 1^ 

^ 2 1 

13 

35 


1^ 

(_5 — 
601 

11 

2101 


) 

l5 

j. 

11 , 

^ T05 ^ 

.4 

^ 2 T 0 


1 

) 


'■ 101 

2101 

X 

+ 

1051 


l2 + _9_ 11 N .3 _R _ J3_ _ 2 _n 

^21 701 35 ^ 101 4201 55^ 

B = 

T? + _!1 + ^ 'l n 4 / _E _1 1 N 

HOI 4201 420 ‘ 601 84 14016 


12 / R_ J3 , _9_ ^ .3 / _E 3 , 15 >) 

^ 21 “ 70 701 ^ ^"101 ” 70 4201 

13 ( ^ ^ Il_ „ _13_ ) ^4 ^ _R ^ _J_ _ _l_v 

^ 101 420 4201 ^ ^ 601 210 1401^ 


t2 ( R_^J3_ _li_ -.5 / __E , ^ 11 S 

^ ^"21^351” 35 6 ^ 101 35 2101^ 

i3 f ^ . 21_ _ 11 ) n4 ( _ 1 _ , _L_ ^ 

101 210 “ 2101 ^ -i- V 210 1051 ^ 


b ) STIPFIIESS CORRECTION MA^TRIX 

0 0 0 0 
0 0 ^ 0 

0 0 M ^ 0 


0 


0 


0 


0 



In this case, we have, 

I 

F = 


M J(R. + 1)^ + Ug^Xl ^ 


F 


M 


(R + 1) -H- 




y 


M 


xy 


M U 

8 

0 


The external force vector Xg is given by 

2 


M (R + 1) XL 


X^=<MU8 


J%2 


> 


0 




(3.14) 


Since the external force Xg is acting at the node itself, 
no transformation is necessary to get the vector Xg in terms 
of nodal coordinates. 

Separating Xg in two parts, we get 


X, 


where 


Q> 


Qg Qg' 


a linear function of displacements U 




a constant vector 


Substituting the values of Xg in equation (3,9), it becomes 
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whGi?e K is the correction stiffness matrix 
Simplifying 

M U + (K - Kg) U = Q 

♦ CjCj 

or M U + K<p U = 0 ' 

^ “2 

K is Called the corrected stiffness matrix: for the case 
2 

being discussed. Elements of the matrix Kg are shown in 
table 3b. 


3,3,3 Expression for Kg 

Figure 4-a shoals the configuration of the forces 
for the case of combined transverse - torsional mode. In 
this dase 5 the body force due to centrifugal action is 
parallel to the Z axis, and the transverse component of this 
force is always zero. Hence, for this case, 


dF = d F = m r dx 

X 

dF = 0 

y 






M (R + x) ^ 
0 
0 


(3.15) 


Substituting this value of Zg in (3,9), we get the equation 
of the type M U + Kg U = Zg' , as before 



TABEE - 4 

a) STIFPliESS CORHECTIOIT MATRIX 


+ J1 ) ^3/ R __3v ^2 . R n 

21 'JO ^ ^ MOi vn-' 1- (oT + ~ 


101 70' 


^21 ■ 35^ ^ I - ^0^ 


1^ (- •r3r:r + 


R . 1 


.) i4( _1_ ) ^3 4 . 

^ 210 ^ Hm 1 V- 


R 


101 105' " 210 ' •“ "101 '420' * 601 


(“ -^T -4J) i2(_£_^) f. R_ -.3 y R 

21 70^ ^ 1 jqJ 1 l-p-, --^j 1 ( — 


21 35 ' 


101 


+ 


,3 y E 


iM — + 


11 

101 " 420 


) (—R + -1 ) i3 / _E 23 ■, 

^601 210 -^ ^ ^ 101 - 210 ' 


4 / I 


1" (: 


210 


b) STIPPRESS COREECTIOU MATRIX K 


0 


0 


0 


O' 


0 


0 


0 


ROTE: Matrices axe related to coordinates, U^, U^. U„. U 

7 57 57 
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whore 


Value of is given in the table 4a. 
3,3.4 Expression for 


Fig. S.4-h shows the force F acting at the end 
B’, Since the forces are parallel to X axis, F = 0 and 


X 


4 


= ' 54 '^ <34 


0 Q!; 

^4 


^4 u + q; 


Here is a null matrix. So when siihstituted in equation 
(3,9) it does not affect K and is given in table 4b. 


3,4 Dynamical Matrix Of The Packetted Group Of 
Turbine Blades 

The packetted group of blades is an assemblage 
of individual blades ajid an overhung band element, connec- 
ted at the free ends of the cantilevor blades. The blades 
and the band element are divided into several segments. 
Nodal co-ordinate system is defined depending upon the mode 
of vibration. The stiffness and mass matrices for an ele- 
ment have been obtained in Section 3,2, Correction stiff- 
ness matrices have also boon derived in Section 3.3. Mass 
and corrected stiffness matrices M and K for the sequence 




NUin6E« OF ELEVEN 


m 'A f % 


lift 4 l\ 


2 w yi-t- 6 


'imLn) f 1 


^ 22143 . 



+ if . 


6tY\4l 

- - - -rirr-*" 

':j4v r 

i, m irt k S 

1 ^ 
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of coordinatos system are obtained -by supor posing and 
Kj^ of the elorrients subject to certain mode of vibration.' 

3.4.1 Tangential Mode of Vibration. 

Figure 3,5 shows a packotted group of n blades. 
Each blade and band aro divided in m and 2 elements res- 
pectively. An additional node is defined at the end of 
each overhanging band. The total number of nodes is thus 
given as 

m . n + (n + l). (3.16) 

It has boen presented that for tangential mode two co- 
ordinates are necessary for each node. Hence total number 
of co-ordinates for the system is, 

2 m.n + 2 (n + 1) . (3.17) 

Tables 5 and 5a show the mass and the stiffness matrices 
for the blade an.d the band elements respectively, 
which arc obtained from the tables 1 ^nd 2 by deleting those 
rows and columns for which the coordinatos are zero. Having 
found and for all elements, superposition is done in 
the following manner. For the first element of the blade, 
first two columns and rows are zero because of the boundary 
conditions. For other elements, the superposition is done 
by adding A of M. (or K. ) to D of M (or K._-,) similarly 
D of (or K^) is added to A of (or 



5^ 


Here A, B, C, D are the sub-matrices of iL (or K^) as shorn 
in the Table 5, This superposition is carried out for 
all the olements of the blade. The ba.nd has 2 n elements 
and (2n + 1) nodes of which alternate node is common to 
the node of the blade, Hote that this corresponds to the 
tip of the blade . 


For each element of the band, (or K^) has 

already been calculated. These are superposed so that A 

and B of M. (or K . ) are added in the following manner, 

1 1 

a) for the starting node of the band, add A corres- 
ponding to 2ran+l, 2mn+2 rows and column 

b) add D of first element of band to D of l‘!^(or K^) 

c) for any element (mn + k) of the band, two cases 
arise 


l) k is odd, then add D of li (or K . ) to 
. mn+k mn+k 


) and A of M, 


D of M , (or K , 


mn+k 


2) k is even, then add A of M , , (or K ) to 
V mn+k mn+k 


D of M , _ (or K , 




) and D of M 
(or K 


mn+k+1 


'mn+k 

). 


Note no change in values of B’s and C*s of I'L (or K^) of 
blade and band takes place during superposition to get M 


and K 
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TiJffGENTIAL MODE 
BMDE PROPERTIES 


STIPEHESS MATRIX 




12EI 

6EI 

-12EI 

6EI 

z 

z 

z 

z 

,3 

P 

3 

2 

1 

1 

r 

1 

6EI 

4EI 

-6EI 

2EI 

z 

z 

z 

z 


1 

1^ 

1 

-12EI 

-6EI 

12EI 

-6EI^ 

z 

z 

z 


n3 " 

2 

3 


1 

1 

1-^ 

l2 

6EI 

2EI 

-6EI 

4EI 

z 

z 

z 

z 


1 

lEERTIA 

1^ 

MTRH 

1 

J1 

111 

9 

-131 

35 

210 

70 

420 

111 

1^ 

131 

-1^ 



210 

105 

420 

140 

9 

131 

13 

-1 11 

70 

420 

25 

1 1 JU 

210 

-131 

-l2 

-111 

1^ 

420 

140 

210 

105 

Note; 


The matrices are j 

lartitioned into four 

(2 X 


sub-matTices , which facilitates the superposition process. 
The partitioned matriz is sho-vm below as, 

A B 

C D 
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TAlGENTIAi MODE 
BAUD PEOPEETIES 

STIEEHESS W-AT-RU. 
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3,4,2 Coupled Transverse - Torsional Mode 

Figure 3,6 shows tho co-ordj.nate system for the 
packetted turbine blades when vibrating in the combined 
transverse torsional mode. The number of node points is 
still given by equa.tion 3, IS, The number of degrees of 
freedom in this case is given as, 

3 m,n + 3 (n + l) (3,18) 

Table 6a and 6b show the elemental stiffness and 
mass matrices for both the blade and the band which are 
derived from table 1 & 2, by properly deleting the columns 
and rows for which the co-ordinates are zero. 

Superposition of these matrices to give 

M & K for the whole system, follows tho same procedure 
described in the last section. 
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A1 


TiiBUS - . 6a 

COUPLED TEiUISVERSE - TOBSION/J, MODE 
bili'le paapppTiEs 


S tif f ne s s I'' e trix 


12EI 


-6EI 


-12EI, 


-'T 


-6EI 


y 


0 


-6EI 




4EI 


0 

-12EI 


r 


-6EI 


0 


j: 


1 

0 

6EI 


j: 


2EI 


1 

0 


JL 


GJ 

1 


1 


6EI 

Z 

1^ 

0 

12EI 
2 

6EI 


JL 


0 


1 

0 

6EI 


4EI 

2 

1 

0 


-G-J 

1 


0 


1 


J1 

35 

-111 

210 

0 

-S 

70 

131 

420 

0 


-111 

210 

jL 

105 

0 

-131 

420 

140 

0 


Inertia Matrix 
0 


0 

J 

X 

3A 

0 

0 

6 A 


70 


-131 

420 

0 

J1 

35 

111 

250 

0 


131 

420 

-li 

140 

0 

111 

210 

1^ 

105 

0 


0 


J 

X 

6A 

0 


0 

5c 

3A 
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ill 


COUPIiEI) TRAWS'VEESE - TOESIOMiOL MODI 
BAl^ro PROPERTIES 


STIEEKESS MATRIX 


12EI 

6EI 


-12EI 

X 

6EI 

X 

X 

0 

X 

^3 

.2 

3 

2 

1 

1 


1 -^ 

1 

6EI 

4EI 


-6EI 

2EI 

X 

X 

0 

G J 

X 

X 

l2 

1 


1 

0 

0 

I 

0 

0 


-12EI 

X 

-6EI 

X 

0 

12EI 

X 

-6EI 

X 

E 




6EI^ 

2EI 

X 

0 

-6EI 

X 

4EI^ 


1 

<M 

H 

1 



-G J 



0 

0 

Z 

0 

0 




DIERTIA MATRIX 


13 

111 

0 

9 

-131 

35 

210 

70 

420 

111 

1 ^ 

0 

131 

- 1 ^ 

210 

105 

420 

140 



I 
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0 

_z 

0 

0 

3 A 
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131 

0 

-11 

-111 

70 

420 

35 

210 

^- 131 ' 

: : - 1 ^ V 

0 

-111 

l 2 

420 

140 

210 

105 


0 


0 



J 


z 


0 


G J 

z 

1 


0 


0 



0 

0 




J 

JL 

3A 


0 



CHAPTEm 4 


RESULTS MD CONCLUSIONS 

In the last chapter v/e have formulated the 
Dynamical Matrix for the t wo modes of vihrations. Respec- 
tive correction matrices are incorporated to account for the 
effect of rotation of the turbine rotor. In this chapter, 
the eigen value problem formulated in the previous chapter 
is numerically solved for both the modes of vibrations. 

4,1 Order of Dynamical Matrix 

Each blade in the banded group is divided into 
5 and 4 elements for tan.gential and coupled Transverse - 
Torsional modes respectivelj?'. In both the cases the band 
spanning any two blade is divided into two elements. Cal- 
culation of natural frequencies for the banded group is done 
for blades numbering one to six. The order MN of Matrix 
for N number of blades for two modes of vibration is given 
by, (Refer. Fig. 3.1, 3.2). 

a) Tangential Mode : ■ 

1^11 = 12 (N) + 2 for 2 

= 10 for n = 1 

b) Coupled Transverse - Torsional Mode 

im = 15 (N) + 3 for n^ 2 


12 


for n = 1 
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For one blade, the values of natural frequencies 
obtained are shown in Table (7), The table also shows the 
theoretically obtained values of natural frequencies when 
the blade is considered as a continuous system. The error 
involved between the two results is in the range of 0*13ii. 
to 5,75% for fundamental frequencj^ to natural frequency, 

4,2 Specifications of the Banded Group 

The following data is used to calculate the 

natural frequency of the banded group 5 used in the H,P, 

Steam 

Section of the 130 MW Ste t n turbine , manufactured by the 
Heavy Electricals (India) Limited, Bhopal. 


1 , 

Length of the blade 

• 

« 

0,88 inches 

2 . 

Cross Sectional area of the 


2 


blade 

• 

0 

0.325 in 

3. 

Moment of Inertia 


4 


-about zz axis 

« 

0.0137 in 


-about yy axis 

# 

« 

0,0089 in'^ 

4 . 

Width of band D 

« 

« 

1*03 inch 

C 

Thiclmess of band B 

• 

♦ 

0,105 inch 

6 . 

Blade pitch, S 

# ■ 

♦ 

0,76 inch 

7* 

Young's Modulus for both 
blade and band 

• 

30 X 10^ lb/ln‘ 

8 . 

Modulus of reaiding for 
blade and the band 

m 

10 X 106 lb/in‘ 
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9, Specific densit 3 i^ for the ; 0,28 Ih/in^ 

material of hand and the 
hlade. 

10, Overhung of the hand : 0,38 inch 

11, Rotor radius at tho hlade root: 17.5 inch 

12* Rotor speed : 3000 rpm 


4,3 Method, for calculation of Eigen Values and Eigen 
Vectors 

MatriJC iteration method is used for the solution 
of the eigen value problem (3,7), The iteration is started 
hy a trial mode^l, 1, ,,.,. 1 "^ ^^rhich is premultiplied hy 
D, The resulting column matrjn: is then normalized hy re- 
ducing the largest numerical element to unity i.e, for 




(4.1) 


where 1 ' and is normalized 

for all i. 


i = 20 is sufficient for the convergence of the normalized 

vector. The normalizing quantity gives the largest eigen 

2 

value for D, which is proportional to 1/p . This gives 
the corresponding minimum frequency, p. 

For obtaining the next dominant eigen value of 

( 28) 

D, Wielandt's Deflation Procedure is used which redu- 
ces the matrix order hy one and eliminates eigen value 
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Corresponding to The iterative procedure gives 

the next dominant eigen value, vrhich in turn gives next 
lowest natural frecraency. Repetition of this method 
gives all the eigen values 3ii order. Rote ej.gen vectors 
are not obtained by this method. 

Since only a few eigen vectors are desired ano- 
ther iterative procedure is used wherein the order of 

(23) 

matrix is not reduced. Sweeping procedure is used to 
supress the unwanted eigen values and the corresponding 
eigen vectors. 

4.4 Numerical Computation 

Computer programs are made for the evaluation of 
natural frequencies for tangential mode and coupled trans- 
verse - torsional mode based on the formulation given in 
Section 4,3, Fortran IV language is used for the programs 
which are run on IBM 7044 - 1401 Digital Computer. The 
details of the Computer Program are given in Appendix D. 

4.5 Results Obtained 

First six natural frequenc-ies for both the modes 
are calculatod for the packetted turbine blades, N= 2 to 6, 
for blade data given in Section 4.2 and are tabulated in 
Appendhc B for f = 0.005 to 0/006. In Appendix A data 
is tabulated for effect of centrifugal force of rotating 
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packetted blades for N = 2 sxid 6 and ^ = 0.02 (0.0231) 
for both tliG modes of vibrations and aro plotted in graphs 
4»1 (a) and 4,1 (b) , For different banded groups frequen- 
cies for fundamental to sixth modes are plotted, a) against 
I'l for parameter 'f and b) against'^ for parameter N. These 
are given in graphs 4.4,1 to 4,4'. 6 and44.3.1 tO/4.3.67'for 
tangential and coupled transverse ~ torsional vibration 
modes respectively. 

4,6 Discussion of the Results 

Mode shapes for first two natural frequencies for 
N = 2 to 6 aro plotted in graphs 4,-2.1 to 4,2,2, 

The results obtained for effect of single blade 
show that -6 l has no effect on the natural frequencies. This 
is quite obvious in the case of very short blades. Results 
for the banded blades, rotor graphs 4,1 show clearly no 
effect of Jl._ on natural frequencies implying the p's for 
banding of short blades are not affected by centrifugal ac- 
tion, Since cases of I, P. and L.P, blades are not considered 
no catogorial statement on effect ofjfLon p'-s is given but 
it is expected that for I, P. and L,P. there should bo pro- 
minent effect. 

Results show that there is marked fall in natural 
frequencies p of banded blades from unhanded single blade. 



Th0 graphs to 4.4,6 are for tangential mode of vibra 

tion of banded blades 5 for bending lias no ef fect whi.ch 
is clear from the node shapes, refer, graph 4,4.1 as mode 
of K .y2 j.s repetition of modes for If = 2 where as for 
p^, i = 2, 3 t ho re is a decrea-sing effect due to If 

for all . Those graphs also show th3,t has marginally 
increasing effect on all p,^, 1 = 2, 3, 4 for banded blades 
If = 2, 3, 4 ... 

The graphs 4,3,1 to 4,3,6 show the results for 
coupled transverse - torsional vibrations. 

Graph 4,3,1 shows that the fundamental frequency 
increases by very small amount as the number of blades is 
increased from N = 2 to 6. This is clear from Fig. 4,4,2, 
for mode shape, which shows that by adding more number of 
blades it puts constraints on the horizontal deflection of 
the overhang3.ng band tips. The graphs 4.3,2 to 4,3,6 show 
that for p^, i ~ 2, 3, 6, there is a decreasing effect 

for the banded blades N =2, 3, 6. These graphs also 

show that has marginally decreasing effect on pj_, 

1 •“ 2 , o , . . , • 

From the graphs of p: against wo note that the 
range of frequency change is more for the case of changing 
the number of blades N than that of . 
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I'fnen the differences in natural frequencies 
i^ 2 are plotted against N, we find that the effect of 
increasing number of blades is prominant in the beginning. 
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TiBLE : 7 

EFFECT OF NUTBER OF BLADES IH TBE GROUP 

Here for an arbitrary ^ and any arbitrary order of 
frequency, the reduction in p is. found for every additional 
balde that is grouped. This is carried out in both the 
types of mode. These reductions for arbitrary ^ = 0.025 and 
p = p3 are tabulated below: 


N 

Cangential 

5 Mode 

P ran s verse - 
Porsional Mode 

2 

313504 

60935 

3 

53001 

5066 

4 

11280 

4269 

5 

67795 

2174 

6 

5368 

1053 
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CHAPTER 5 


COI'ICLUSION Ai'ID REGOlCiMfDATIONS 

In the present work, the problem of banded blades 
is analysed for its vibration, tangential and transverse - 
torsional modes* Finite element method is used. 

The results for a short blade show that centrifu- 
gal forces have neglegible effect whereas banding have 
appreciable effect. In addition to method of computing 
natural frequencies and mode shapes of banded blades, data 
has been calculated and presented which can bo useful for 
the design of the turbine blade of H.P. Section, The rela- 
tive merits of the method can not be- presented due to the 
fact that only a limited work has been published for the 
vibrations of banded blades_.in particular the example taken 
in the earlier work is different from that of the present 

■■ t 

work. 

It is recommended that hinged flexibility and 
stagger angle at the roots of the banded blades of H.P. 
Section bo incorporated and results compared x^ith earlier 
vrork. TW.s vr.lll give a clear idea of the relative merits 
of tho two approaches. In addition to this, problem of 
banded blades of I.P. Section can be solved wherein twist 
of the blade can be incorporated* Since L.P, Blades are 
connected by lashing these can also be solved by the above 
method incorporating, in effect nf ; centrifugal forces as well 
as taper twist* .x' - . 
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APPEiroix: A 


EEEECT OF ROTOR SPEED OR THE VIBRATIOF PEEQUEITOIES FOR 

BiRPED BliilES 

In this appendix, the results for fundamental and higher 
order natural frequencies are tabulated for "various rotor speeds. 
Two cases of banded groups having 2 and 6 blades are considered for 
=5 0,02 (corresponding--*? = 0.231 ) for tangential and coupled 
Transverse - Torsional inodes of vibrations. 


TABIE - A.1 

TAKfGENTIiSJ, MODE OF VIBEATIOR 

= 0 . 020 ) 


p 

■^1 ^ 

“I 77“ 

1 15 

T 

f 30 

I « 

T 60 

5, 60 

Pi 

17893.2 

17893.2 

17893.2 

17893.2 

17893.2 

P2 

4 9930,2 

49930,2 

49930.2 

49930,2 

49930.2 

P3 

114221.4 

114221.4 

114221.4 

114221 .4 

114221.4 

P4 

120678.1. 

120678.1 

120678.1 

120678.1 

120678.1 

P5 

187136.5 

187136.5 

187136.5 

187136.5 

187136.5 

Pe 

198738.5 

198738.5 

198738.5 

198738.5 

198738.5 

Pi 

17893.2 

17893.2 

17893.2 

17893.2 

17893.2 

P2 

24225 

24225 

24225 

24225 

24225 

P3 

36777.7 

36777.7 

36777.7 

36777.7 

36777.7 

P4 

49921,8. 

49921 .8 

49921.8 

49921.8 

4 9921*8 

P5 

61201,4 

61201.4 

61201 .4 

61201.4 

61201,4 

% 

68807.4 

68807.4 

68807.4 

68807,4 

68807.4- 


imjE - A. 2 
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OOlIPIiED TRAMS'VEESE - TOESIOffii MODE 

(f = 0.0231) 



Pi 

20472.6 

20472.6 

20472.6 

20472.6 

20472.6 

P2 

24058.7 

24058.7 

24058.7 

.24058.7 

24058.7 

P3 

45136.7 

45136.7 

45136.7 

45136.7 

45136.7 

P4 

66982.6 

66982.6 

6 6 982.6 

66 982.6 

66982.6 

P5 

99365.4 

99365.4 • 

99365.4 

99365.4 

99365.4 

P6 

133408.4 

133408.4 

133408.4 

133408.4 

133408.4 

Pi 

20540.3 

20540.3 

20540.3 

20540.3 

20540.3 


Pg 

20963.8 

20963.8 

20963.8 

20963.8 

20963.8 

P3 

22895 

22896 

22895 

22895 

22895 

P4 

26638,3 

26638.3 

26638.3 

26638.3 

26638.5 

P5 

33185.7 

33185.7 

33185'.7 

33185.7 

53185.7 

1 

Pg 

43015.5 

48015.5 

48015.5 

48015.5 

48015.5 - I 

■ ", I 



IPPEEDIX B 


DATA FOR NATURAL FREQUENCIES OF VIBEiTION OF BINDER BLADES 


In this eippendix values of natural frequencies for different 
modes for a stationarj?' banded group of blades in both tangential and 
coupled transverse - torsional mode is given. The number of blades N 
is varied from 2 to 6. For every group 5 the band to blade stiffness 
ratio ^ is varied from 0,005 to 0.06. The values of natural frequen- 
cies p's are tabulated below. 


The data is presented in three sections. Sections B,1 and B,2 
give results for tangential and coupled transverse - torsional modes 
respectively and in Section B.5 coordinates for first two eigen vectors 


for both modes for N = 2 to 6 are given. In each sections B, 1 and B,2 
six tables are given for first six natural frequencies. The corres- 
ponding natural frequency p for the single blade is given on the top 


of idle table, 


TABLE B1.1 

FUND/lEmL NATUEiJb FEEQUENCy 

(F . V ^ = 24,304.34 c.p.s) 

Single blade 



— ^ f 
^ 1 

“T — ^ 

■"■T 

5 

T 

i 

6 

0.005 

16879.25 

16879.25 

16879.25 

16879,25 


16879.25 

0.010 

17238.36 

17238.36 

17238.36 

17238.36 


17238-36 

0.015 

17575.61 

17575.61 

17575.61 

17575.61 


17575.61 

0.020 

17893.25 

17893.25 

17893.25 

17893.25 


17893.25 

0.025 

18193.10 

18193.10 

18193.10 

18193.10 


18193.10 

0.030 

18476.72 

81476.72 

18476.72 

18476.72 


18476.72 

0,035 

18745.49 

18745.49 

18745.49 

18745.49 


18745.49 

0,040 

19000.60 

19000,60 

19000.60 

19000.60 


19000.60 

0.045 

19243.13 

19243.13 

19243.13 

19243,13 


19243.13 

0.050 

19474.03 

19474.03 

19474,03 

19474.03 


19474.03 



miEB.1.2 

SECOID higher NITUE/lL erequeetcy 
’ single blade ~ 152,39.30) 


1^" 

1 — 

5 

1 — 

i 5 

T— T 

^ ^ 

T 

1 

6 

0.005 

49781.5 

36404.0 

2 96 93.0 

25882.0 


23530 

0.010 

49836.2 

36539.7 

29875.0 

26100,8 


23770 

0.015 

49880.3 

36662.5 

30046.7 

26307.7 


24003 

0.020 

49930.2 

36779.0 

30209.9 

26504.3 


24225 

0.025 

4 9 977.2 

36889.8 

30365.3 

26691.3 


24435 

0,030 

50021.6 

36995.4 

30513.4 

26869.5 



0.035 

50063.8 

37096.3 

30654.9 

27039.6 



0.040 

50104.0 

37192.6 

30790.2 

27202.1 



0.045 

50142,3 

37284.8 

30919.6 

27357.5 



0.050 

50178.8 

37373.0 

31043.6 

27506.2 





T/iEIE B.1, 

,3 





THIRD HIGHER HITDRAIi EEBQUEI3CY 





(p 

■’^single blade 

427,725.9) 




V 

/ 

1 ^ 

IT ^ 
i ^ 

^ A 


T 

i. 

6 

0,005 

90568.1 

61198,4 

49772.7 

41848.6 


36385.4 

0.010 

110196.7 

61201.8 

49841.3 

41968.74 


36557.9 

0.015 

11306 9,0 

61211.4 

49892.5 

42029.5 


36661,6 

0.020 

114221.4 

61220.3 

49940.6 

42145.36 


36777.7 

0,025 

115021.4 

61228.9 

49986,1 

42226.8 


36878,0 

0,030 

115687.2 

61236.9 

50029.2 

42304.3 



0.035 

116283-.4 

61244.6 

50070.2 

42378 



0.040 

116836.3 

61251.8 

50123.7 

42435.8 



0.045 

117358.2 

61258.3 

50161.05 

42500.8 



0.050 

117856.2 

61264.8 

50196.6 

42564 .*9 




POUBIH HIGHER HATURjili FEEQDEIKY 
blade = 845584.6) 


1 

2 

I 

3 

~T 1 

1 ♦ 1 

5 

T 

i_ 

6 

0,005 

93529.6 


90412.6 

65532.0 

57009.3 


49794.7 

0.010 

116757.6 


109724.8 

65539.4 

57006.2 


49817.21 

0.015 

119676.0 


112522.6 

65521.1 

57024.4 


49886.1 

0.020 

120678.1 


113656.1 

65542.4 

57057.5 


49929.8 

0.025 

121664.0 


114443.4 

65154.3 

57077.0 


49970.4 

0.030 

122173.1 


115100,3 

65544.6 

57082.2 



0,035 

122762.4 


115688.1 

65545.1 

57098.9 



0.040 

123322.3 


116232.7 

65527.0 

57135.7 



0,045 

123854.3 


116748.0 

65527.0 

57151.3 



0.050 

124358.7 


117238.6 

65528.3 

57166.9 






TABIF, B. 

1.5 





FIFTH HIGHER HATURAE FEBQTIEIKY 





Y 

single blade 

= 1,403,307) 




H j 

, '* , 

! ^ 


3 

_1 j 

1 * } 

5 

T 

6 

0.005 

115780.6 


99523.1 

90439.1 

67581.34 


61170.8 

0.010 

137410.1 


113S16.0 

109766.0 

67610.94 


61185.5 

0.015 

101132,1 


116823.0 

112604.0 

67619.9 


61195.5 

0.020 

187136.5 


118007.0 

113856.7 

67683.3 


61201.4 

0.025 

192635.4 


118838.0 

114600.9 

67602.4 


61204.3 

0.030 

195661.8 


119561.6 

115194.0 

67616.4 



0,035 

195246.2 


120190,6 

115808.6 

67619.4 



0.040 

195757.9 


120769.6 

116823.5 

67618,5 



0.045 

196 2 5 3.3 


121314.0 

116801,8 

67618.8 



0.050 

196556.0 


121821.5 

117481.7 

67618.0 
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B'l.e 

Sn:iE HIGHER NATURAL PHEQUEHCY 
^^single blade 2,333,891 ) 


N 1 

I 

2 

i 

3 

1- ^ 

H 

5 

T 

6 

0.005 

146532.0 


101088.4 

94487.9 


90088.8 


68113.5 

0.010 

183257.3 


124805.5 

113797.4 


108770.1 


68185.0 

0.015 

193371.7 


105618.2 

116219,6 


111281.9 


68782,7 

0.020 

198738.5 


126493.2 

117729.4 


111957.7 


68003.4 

0*025 

211326.5 


127153.0 

118468.5 


113389.4 


68807.4 

0.030 

213876.4 


127658.0 

119392.0 


113648.2 


69487.3 

0.035 

214169.2 


128182.6 

124392.1 


114127.3 



o.ow 

215561.6 


128702.1 

120059.7 


114437.5 



0.045 

210359.9 


129176.8 

121.120*5 


114537.0 



0.050 

216638.4 


129654.5 







» t- • I 


I'UNDAEvEM'M, HATDEAX FEEQUEHCY 
' '’stogie blade = 50184.29) 



TiiBIE B-.2.2 

SECOHP HIGHEE I/sTUBil EEEQTJEKY 

(p . , , , ^ = 36793.8) 

single blade 




A 


0,0077 

58068.3 

40696.1 

30719.8 

26230.2 

24052.0 

0.0231 

45136.7 

32182.2 

26687.9 

24242,6 

22895.0 

0,0385 

40553.5 

30070.2 

25801.8 

23827.5 

22674.85 

0.0539 

38184.7 

29076.6 

25391.6 

23630.6 

22458.3 

0.06 91 

36713.1 

28465.2 

25127.1 

23493.1 

22281.5 


TABEE B.2.4 

BOUEffl HIGHER HATTJRili EEEQUEHGY 

(p . ^ ^ = 189375.9) 

■^single blade 


BBI 


■ f , 

1 ^ 

1 

1 


0,0077 

74951.2 

61818.0 

51194.3 

40380.5 

33344.61 

0^0231 

66 982.6 

48814.4 

38378.25 

31611.26 

26638.3 

0.0385 

6 0954.8 

43568.7 

34674.14 

29449.3 

25453.79 

0.0539 

56574.4 

40751.6 

32869.40 

28441.69 

24723.0 

0.0691 

53326,3 

38 975.39 

31767.11 

27828.6 

24025.7 
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B.2.5 

FIFHI HIGHER HATURAI EEEQUENCY 
^^single blade~ 210847.0) 



T 

i 

2 

1 — r 

-I ^ 

X 

± 

4 

T 

1 

c H 

^ 1 

6 

0.0077 


113539.6 

73071.6 


64848.8 




56802.3 

47 922.4 

0,0231 


' S9365.4 

66917.5 


50700.6 


42554.8 

33180.7 

0.0385 


94108.3 

62109.4 


45162.9 


37885.09 

30825.3 

0.0539 


90439.2 

58549.19 


42122.4 


35537.9 

29507.2 

0,0691 


86805.3 

55746.0 


40187.2 


34091.7 

27953.4 


TABIE B.2.6 

SKTH HIGHER NATURAL EEEQTENCY 
(Psii-gle blaae = 


1 

2 

1 

5 V 

T 

1. 

5 

X 

JL 

6 

0,0071 


141927.7 

79056.8 

71051.1 


65181.3 


59839.8 

0.0231 


133408.4 

67621.5 

66762.58 


51772.0 


48015.5 

0.0385 


105666.4 

60037.49 

59062.13 


96099.5 


43672.3 

0.0539 


97314.15 

54960.8 

53795.64 


42941.15 


40012.7 

0.0691 


90988.05 

51243.5 

50155.34 


40917.6 


37833.25 


T/iBIE B.3 
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MOIE SHIPES 

(a) First two eigen ■vectors calculated for tangential inodes of vibrations 
of banded blades 1 = 2 to 6 are given below for = 0.025 
a) ]\| = 2 



\ 0.041, 

0.447, 

0.149, 

0.758,' 

0.300, 

0.937, 

0.472, 

1.000, 

0.647, i 

^1 - 

0.041, 

0.447, 

0.149, 

0.758, 

0.300, 

0.937, 

0.472, 

1.000, 

0.647, : 


0.662 , 

-0.493, 

0.662, 

-0.493, 

0.662, 

-0.493 




i 

1 0.067, 

0.698, 

0.222, 

1.000, 

0.400, 

0.973, 

0.554, 

0.766, 

0.672, : 

III 

-0.067, 

-0.698, 

-0.222, 

-1.000, 

-0.400, 

-0.973, 

-0.554, 

-0.766, 

-0.672, 4 


.803, 

-0.348, 

.000, 

-0.000, 

-0.803, 

0.348 ( 

r 


1 


b) 

N = 3 









0.041, 

0.448, 

0.150, 

0.758, 

0 . 301 , 

0.938, 

0.473, 

1.000, 

0.647, ! 


, 0.041, 

0.448, 

0.150, 

0.758, 

0.301, 

0.938, 

0,473, 

1.000, 

0.647, 1 

IT = 

1 - 

0.041, 

0.448, 

0.150, 

0.758, 

0.301, 

0.938, 

0.473, 

1.000, 

0.647, i 


0,662, 

-0,494, 

0.662, 

-0.494, 

0.662, 

-0.494, 

0.067, 

0.699, 

0.222, ? 

1 

1 0.054, 

0.577, 

0.186, 

,0.901, 

0.358, 

1.000, 

0.531, 

0.947, 

0.689, C 


kw 

0.000, 

0.000, 

0.000, 

0.000, 

0.000, 

0.000, 

0.000, 

0.000, 

0.000, ( 

V ss£ 

2 ~ 

-0.054, 

-0.577, 

-0.188, 

-0.901 , 

-0.358, 

1 

• 

O 

o 

o 

-0.531, 

-0.947, 

-0.689, -( 


0.758 , 

-0.494, 

0.379, 

-0.231, 

-0.379, 

0.231, 

-0.758, 

0.462 
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1 = 4 









= 1^5 A, A, A, B, Bs 

B, B, B 

7 

V 







Ytoere A =,| 0,041, 

0.447, 

0.149, 

0.758, 

■ 0 . 300 , 

0.937, 

0.472, 

1.000,; 


0 . 672 , 

0.969 








b' =1,0.662, 

V 

-0.093 








= ^A, B, -B, -A,’ 

0.735, - 

-0,492, 

0.519, 

-0.348, 

0.000, 

0.000, 

-0.519,; 


0.348, 

-0.7351, 

+ 0.492 







Where A = (0.049, 

0.681, 

0.532, 

0.175, 

0.059, 

0.3419, 

0.999, 

0.519, 

1.000, 


0.935 \ 








B =1 0.026, 

0.220, 

0 , 072 , 

0.355, 

0.141, 

0.414, 

0.214, 

0.414, ! 


0,285, 

0.5873 1 

- 







N = 5 










By By 

B, B \ 

J 





I 


Where A = as detailed for 

in tlie case 

N = 4 



: 1 


B = as detailed for in the case 

1 = 4 




"2 

“( ® 

1 

) 








^Hhere A 0*0466, 

0 * 500 , 

0,165, 

0.820, 

0 . 325 , 

0.976, 

0 . 501 , 

1.000, 


0.673, 

0.948 

i 







B =f 0.028, 

0.309, 

0.1029, 

0 . 507 , 

0.201, 

0.605, 

0 . 309 , 

0.618, 


V ■ 

0.416, 

0.5859 

r 







0 = ^0,0, 0 

it 0 y 0 • 0 y 

0 . 0 , 

0.0, 0 

'.0, 0.0, 

0.0-, 

0.0, 0, 

.0 ] 



N = 6 


j -“s A, A, A, 3 ^ 3 ^ ^ ^ 


^ j A j 

?/he re A 
B 


defined for 7^ in ilie case If = 4 
as defined for 7^ in the case IT = 4 


7 , 


./;., E, U, - 0 , -B, -A, 0.069, -0.493, 0.600, -0.427 , 0.346, -0.246, 

0.0, 0.0, -0.346, 0.246, -0.500, 0.427, -0.693, 0.494 I 


B 


III 

o 

• 

VJI 

0.483. 

( 

m 

o 

o 

0.665, 

= j 0.032, 

0.353, 

0.732, 

0.486, 

= .1 0.012, 

0,129, 

i 

0.267, 

0.178, 


.232, 0.705, 0.360 


C 
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B ~ 0,00, 0,00, 0,000 — 0,00 

If K 4 

A, A, A, A, -0,618, 1,000,-0.09, -0,601, 1.000 


-0,001, 

-0.601, 

0,999 

-0.618, 

1.000, 

0.090 

Where A = 

-0.053, 

0.463 


-0.381, 

0.930 

1! 

, -A, 

—0,652 


-0,703, 

0,000 


-0,703, 

0.652 


N = 5 


•0.000, -0.601, 

1.000, 0.001 

0.000, -0.192, 

0.768, 0.000 

0.000, -0.593, 

0.976, 0.000 

1.000, -0.167, 

-0.438, 0,707 

0.000, -0.991, 

0.438,-0.707 

•1.000, -0.167 



1 

— A, A, A, A, —0.618, 

1,000, -0.090, 

-0.601, 0,999 


-0.001,-0,600, 

0.999, -0.000, 

-0.600, 0.999 


0.000, -0.601, 

0.999, 0.001, 

-0.618, 1.000 


, -0,090 




Whem A =-0.053 , 0,463, 

-0.000,-0.192, 

0,768,-0*000 


-0,381, 0.930, 

-0,000, -0,593, 

0,976,-0.000 

2 

= A, B, C, -B, -A, -0.640, 

-1,000, -0.139, 

-0,496, 0,800 


-0.456, -0.189, 

0.308, -0,734, 

0.189,-0,309 


-0.734, 0,496, 

-0.809, -0.456, 

-0.640,-1,000 


-0*139 
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(b) First two eigen vectors calculated for coupled transverse 
torsional modes of banded blades N = 2 to 6 are given 
below for = 0,0385 

N = 2 


7^ = A, A, -0,618, 1.00, -0.093, 0.601, 1.00, 0,00, -0.618 

1.00, 0.09 

Where A = . 0.053, 0,046, -0.000, -0,192, 0.768, 0,000 

-0.381, 0,930, -0.000,-0.598, 0.976, 0.000 



-B, 0.377,-0.427, 0.322, 
0.427, 0.322 


-0.000, 0,000, 1.000,-0,377 



Where 

B = 

0.028 

,-0.239, 

0.060, 0.096,-0,367, 

0.111 




0.183 

,-0.403, 

0.145, 0.271,-0,395, 

0.157 


N = 

3 





V 

1 

11 

s J 

-0.618, 

1.000, 

-0.09, -0.601, 1.000, 

—0,000 




-0.601, 

1.000, 

0.000,-0.618, 1.000, 

-0.09 


Where 

A = 

-0.053, 

0 .463 , 

0.000,-0.192, 0.768, 

0.000 




-0.381, 

0.930, 

-0.000,-0.539, 0.967, 

- 0.000 

^2 

= A, B 

j -A, 

0,556, 

-0.819, 

0,187, 0.260,-0.410, 

1.000 




-0.260, 

0.410, 

1.000,-0.556, 0,819, 

0.187 


Whe re 

A = 

0.049, 

-0,422, 

-0.000, 0.173,-0.669, 

-0.000 




0.333, 

-0.766, 

-0.001, 0.503,-0.774, 

-0.001 



appendix g 


COMPUTER PROGRAM 


This appendix gives the description of the computer 
prograitis used to calculate the natural freQuencies of the 
handed group of blades. Two different programs, TAi'IGSN and 
TRATOR are made for tangential and coupled Transverse - 
Torsional vibration modes. A listing of the programs and the 
subroutines used are added at the end, 

C.l Main Program TA'I GEN 

This program calculates the natural frequencies of 
the tangential mode of vibration for the packetted steam tur- 
bine blades. The data required is to be defined in the program 
itself. Prints out the natural frequencies and the mode shapes. 
Only the natural frequencies are calculated if MODE = 0 is 
defined. Subroutinos used are given below. 


0,1,1 Nail 

.0 of Subroutine 


STIFF 

a) 

Argiiment list 

m 

A, D, M 


A 

" ■ 

Length of the element 


D 

W' 

Stiffness Matrix 


M 

s 

Number of element 

b) 

Common List 

• 

« 

ZBl, E, R, W22 


ZBl 

t 

MJ- of the section about 


E 

1 

axis of the blade. 

Youngs Modulus 
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• Distance of the inner end A 
from the rotor axis. 

¥22 : This is square of the rotor 

speed. 

c) Dimension Statement 

D ! (14, 14) 

d.) Purpose : To calculate the basic stiff- 

ness matrix of the matrix 


C,l,2 Kamo of the Subroutine:' MASS 


a) Argument list 

A, M 
C 

b) Common Statement 

RHO 

H 

AB 


A, C, M 

Def jjied as earlier 

Mass Matrix of the element 

RHO , H, AB 

Mass per cubic inch = 

ZT/g for blade 

Elemental length 

Area of the blade cross section. 


c) Dimension Statement 

C : (14, 14) 


d) Purpose 


To establish the basic inertia 
matrix used in TAI'fGEN 


0,1.3 Name of the Subroutine: MATBI . 

Argument List : A, N . 

As Matrix to be inverted, and the 
inverted matrix itself 

N s Order of A 
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b) Common List : 

c) Dimension A : 

<i) Purpose : 

C.1.4 Name of Subroutine : 

a) Argument list ; 

E, RHO 

I, J, K, L : 

b) Common List : 

AM : 

AK : 

c) Dimension AM : 

ak: : 

d) Purpose ; 

Ctl*5 Namo of the Subroutine i 
a) Argument List 

AK, AM : 
mi i 
MO } 

X i 


NIL 
(N, N) 

To Invert the matrix, 

BAND 

E, REO, I, J, K, L 

Youngs Module & specific 
mass density for the band* 

Co-ordinates of the band, 

AJ4, Mi 

Mass Matrix of the whole 
system 

Stiffness matrix of the whole 
system 

The value given for MM. by 
equation 3,17 

Same as MM 

To couple the blade matrices 
with those of the band 

EIGMOD 

AK, AM, Ml^, MO, X, EIG 

Mentioned earlier 

Order of matrix AK, AM 

Number of lowest natural 
frequencies required. 

The matrix of mode shapes 

Natural frequencies 



b) 

c) 
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Common List : 

Dimension Statement: 

M, AM : 

X : 

EIG : 

d) Purpose : 


MIL 

(M, M) 

(W'J, m) 
iVD) 

To calculate the natural 
frequencies and the mode 
shapes. 


C.l.G Name of the Subroutine ; SWEEP 

a) Argument List : AK, AM, D, MN, III, X 

AK, AM, M1-I-. ,! Mentioned earlier 
D : Dynamical matriJC 


III : Number of lo\-7est eigen value 
being calculated 

b) Common List NIL 


c) Dimension Statements 


Ali, AM ; (M, MI) 

D : (MN, MN) 

X : : (MN, MO) 


d) 


Purpose 


: To sweep the dynamical matrix 
by one dimension. This is 
called from EIGMOD, 



Ill 


C*2 He in Pro grain 


TR/ilOR 


This program calomlates the natural frequencies 

Of the coupled Transverse - Tor^lon = T >•' ^ 

iorsional Hodes of the banded 

turbine blades. The date, nr -i-i i . 

aata of the blade specifications is to 

be defined In tho programme Itself. 

STIFF 
A, D 


0.2,1 Name of the Subroutine 
a) Argument List 

A 

D 

t*) Common Statements i 

RliO, H, AB J 
E, R 

XI ; 

I : 
ZB2 j 

, c) Dimension Statement* 

D 

d) Purpose ; 


Elemental length 

Basic Stiffness Matrix 

RHO, H aB, XI; ZB2, B, R 
lif22, I. , ’ ■ ’ ^5 n, 

Mentioned Earlier (for the 
blade element) 

Torsional inertia constant 

Torsional stiffness constant 

M.I of blade section about 
the I-I axis. 


(6j 6) 

To calculate basic stiffness 
matrix 


C.2,2 Fame of the Subroutine : MASS 


a) Argument list 


A, C 


A, C j Mentioned in C,l.2-a 
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b) Conaiion list ; RHO, H, AB, XY 

, H, ABjZY: Mentioned earlier 

c) Dimonsion Statement: 


C : (6.6) 

C,2,3 Name of the Subroutine : BMID 

a) Argument List : E, RHO , I, XY, II, 12, 13, 

Jl, J2, J3 


E, RHO, Y, XY 

■"1, 12, 13, : 

Jl, J2, J3 

b) Common Statement : 

AIC, AM ; 
S : 


Mentioned Earlier 

Coordinate of the band 
element 

AK, Al^, S, ZP2, AP 
Mentioned earlier 
Band elemental length 


ZP2 : M.I. of the band section 
about X-X axis. 

AP : Cross-sectional area of the 
band , 

c) Dimension Statement: 

M i (MN, xMlO 

AM : : i , M'j) 

tl) Purpose s To couple individual blade 

mass and stiffness matrices. 

G,2,4 Hame of the subroutine : MATIIi 
Described in 0,1.4 
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C.2.5 Name of the subroutine 

a) Argument list 

AK, MO, 3IG 

im 

b) Common list 

c) Dimension 

iiZ 
31 G 


EIGEN 

Alf 5 M, MO , EIG 

All mentioned earlier 

Order of matrix in the 
case of combined transverse- 
torsional mode as given in 
the text (Equation 3,18). 

NIL 


(M, M) 

(MO) 

This Calculates only the 
lower I'D natural frequencies 
and no mode shapes. ■ 


d) Purpose 
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$ IB JOB 

SiHFTC TFAfOR 


C 


^ »,' .](. 

C 

C 


’*■' * '“' ■''■ '''■ * ■''' * * * * * * * * * -jf '!{■ -X- -X X -X * -X- -X -Jf * * * * -Jf- * * * -i 

C !->PObiiAM |v!AuE BY, VASAnT J. BHIDE , M . TECH . ( mECHc ) 

rxxxxxxxxxxxxxxxx . ... J : ! 7/8/1972 

^■' R' R'*’ ^ 7C- yi- -A- -)c •?(■ -'v- -yr <- -^C if "X' f<r 3f ^X- -X- -;f -X- -X- 'X* -)(- -x- ^ 

c 
c 
c 
c 


c 

c 

c 


Coupled tpahsvEpse and topsional vibrations of h„ po turbine bl 

D I MLmSIoM Ak ^93?Q3 ) f Am (.93 ,93 ) ,Am 3( 18, 18) ,AKb( 18 ,18 ) ,STlF( A,6 ) 

■X , /vf.lASS ( 6 'J 6 ) , C ( IB ) s E I G ( 6 ) ? XXX ( 93 ? 6 ) 

C 0 / : I C N / A A A / i'? H (,) <1 H 9 f \ B s X Y 
*'■ fjf ifvtC'i! / / / b B 5 E 9 R , iiV 2 2 9 7 

Cor/MON /CCC/AK^AMsSsZPl sAp 

call. ELOV(320Cj) 

CALL Flunk 32wOu) : 

MODE=0 , I 

physical PPOPRETIES OF THE BlADE AnD THE BAnD ^ 


DATA E » ITmAX 9Cj/300u0000o 9 1 8 9 12000000 0 / 

RHO=Oo28/3a6o4 

A=0„88 

P A 'ss 1 7 

AB=212o857/60Ao2 

AP=2oP67*26oi62/65492 

ZP1s:26<,16 2/25„A/12o-»K 2o667/25o4)**3 

ZP2«2 9667/25oA/12o*( 26o162/25o4)**3 

ZB2«5a63aA/65A.2/654«2 

5=19«301/2S»4 

S = S/,?. 

B=13.75/28«4 

D=17^36/23«4 

PPInT l0»Aib»D9AB9ZD29AP9S9Zp] 

Print ii9E»C9»rho 

c^c L 19 U -louivalEnI bpEadth a^d thickness of THE section 

P|=3«141393 

w?=o« 

M = 4 
AflaM 
H®A/AN 
C 

c Topsional Constants eval^Jted 
T = l« 

sum*o« 

114 X=T#PI»D/8/2o 

MZ = 0 

112 MZ*MZ+1 
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IIS’ 

X--=X/2.. 

I F(XoGTo4.> ) Go TO 112 
X--!. XFMX) 

OO 113 , 

>=X-X 

X - (X-i. /X ;/ (X + i., /X) 

X^^X/T«)f5 
6lJM=SUM+X 
T T + 

If- (XoGT,.OJa)uul ) GO To 114 
X :. { 1 c - 1 '■! 2 .. / P 1 -:(-B / D* SUM 1 / 3 o 

Y:-X l-'-G/H 
XY = Xj;!:NiC<.<H/3o 
pr: rr'T i;.:,Y,XY 

ccc 

c InY'PL ..-riFF.. 9 mass matrices calculated ( — STo 110 ) 
c To account EFFECT oF Co E^REMOVE CARD n/cooo 104 
CCC 

Do 100 
Do 10 ( = 

AkU < I » J ) -Wo 

1 uO l»J!“Oo 

W22 = V./2#W2 
DO 110 I = 1-.M 
Al=f-1 
I3=3*n-1) 

R=RO+HKA{ 

lw 4 I F ( 1 T., 1 ) 00 To 102 

CALL ST II F(H, ST IF ) 

CALI, MASS (H? Amass ) 

102 DO no j = n6 
J3= I 3 + J 
Do 1 10 I =] .6 
|.:' = f 3 H. 

AkO( J3»l 3)=:AKB( j3»LS)+STlFl JsL) 

AMBI J 391. 3 ) -AMl'il J3 »L3 ) +AMASS ( J 9L ) 

'CCC ' 

c boUfviDApY CONDITION'S APPLIED ( — 120 ) 

M3*MM-3 . X ■ X VT 

DO 12U I-1,M3 
DO 120 J=1»M3 
AKtM n,J)=AKB( H-3»J + 3) 

I 2 O AMB { I » J ) = AMO (I + 3 » J + 3 ) 

C ■ ' ■ C :T 

DO 2550 M-6»10 

C , ■ ■ X :x’ 

RATIO^OoDAO 

IF{N«.r:o.l} GO To 111 
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28 0 / r- 2 = Z b 2 ).'■ R A T J 0 

D = 8r,RT ( 1 ,’„';2.p?/AP) 
B = Af>/D 

V 2/12 , 


Lai **" i' i 

1114 X " 1 I »-Li/ L /,? „ 

1112 .'■'/■•M/'a 
X“X/2,. 

I M ^ >iCiT 4 ) Gn To 1112 
A = EXrMX) 

00 1113 1=1, MZ 

1113 y-'XnX 
Xt ( X—K /X ) / { X + i ^ /X ) 

X = X / T •'! 

S0M=4Ui''4j-X 
T=T.f.2„ 

if- (Xu T o I u u 1 ) Go T 0 1 1 1 4 

y = ( I ~ I P 2 i. / P [ iHf- 8 *-B / D StJM ) / 3 
X l = X 

Y=xj»r./H 
XY = X J*RHrHfH/3o 
n*Y,XY 

111 continuf 
ccr 

MH=3»M*M+3 

1 F ( 0 1 ) MN*MN+N+N+N 
DO 1 30 !*=I,MN 
Do 13U J=1,M0 

^ AKilfJlaUi, 

13‘-' AM{ I » J ) =1’ a 

M5at,M4LfM0 0^' ,1A5S Af.iC- SFIFmEOS h^TRUCES GOR THE SySTEm. 
DO 140 1*1, |\j " 

N3a3#M» { I-l ) .. 

Do J 40 J» 1 ,'«3 
MJ-J4-N3 
DO 14') K»UM3 
fix 

AK(N,J,MK) ~AKPIJ,K) 
lA 0 AM(WJ*'M)<S«A*4B( J,K) 

etc Coupt.iNCi uu::, lo band ( « o o oo 150 ) , ' 

I f (M»EQ.. I } GO lo 160 , . : . 

! 2«n~i ' ■ . 2, ■ 

DO 150 t I »1,N ■ 2'.. ' ■ ' ' ' ' ' 2 

Il»U+3 ' '2' ^ . 

. 2 R," , 



J'X) 
16 0 
ccc 


18 0 

170 

C 

85 

84 

16 

C 


,,255 

.1 


I 3= T 3 + 3 


r I 
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1 1 » 1 2 9 1 3 s J 1 j 2 5 j 3 ) 


J2 = J'1-1 
.i3 = J2-l 

C A L. I f A i) ( r •) R Mo , Y , X Y 
K1*JI 
rvi-'-jR 
K3 = J'3 
i, 1 = f 1 + 3 
L2-^T;'! + 3 

I. 3--- 1 3 + 3 ■ 

i AM3 n.. 9 r.Ho s Y <»XY 9 K1 9 K2 9 K3 9 L 1 9 12 9 L3 ) 

CONTI NUi: 

dynamical. ,MTpIX For the system oBTAUjED U ST m/ 
IF(N„.Foa ) MiM = W 3 • ^1= N/o- 

CALL MATIN (AK»MM) 

DO 170 I=1,mN 
DO IBO J=1,MN 
C{J)5^0« 

Do 16u L=l9MM 

r { .1 ) =C. ( J ) +AK ( I 9 L ) ^^AM { L 9 J ) 

DO 170 L. = 1»MN 
C I L ) »1000000o0 30«»C ( L ) 

A K n ♦ L ^ ■ C ( L ) 

PRINT IfNfRATlO 
M2«6 

IF(mO) 84»84»85 

CALL *:IGM0D C AK i Am »MN »6 »XXX 9 EIG ) 

GO To 86 

CALL eigen IAK*MN»M2»EIG) 

Conti nue 


■170 


IF{N*CO#I J GO To 2550' 
RATI0«RATI0+0,,U1 

IF(f?ATlO»L£»0«O4?)/-6O TO 280 • 

coNT rwuE. ... „H 

ForwAT' H-%m **N0«0F''.96't AO£i-;';;+ti 


20Xy*8LADE To BAND STiFFo nATIo 


iiiicc:ai 

^r- i L H- ) 



2 

f' r ( .'■ 

.. 

^MODE,, 

3 

; f ! ( lu 

X»^f 

^ 1 . i iw V 

iliilil: 


K* Iv 

/ E 2 « 3 

ID 

i O'.'MA !' ( /.■ 

0 >- . . 

I' A- 

li 

format (1 



L 2 , 

* — 

I ^x- 


„_^||||li|S 


1^ 


14 

(' :■! -’A I ( / 

IW 

'■ ! ‘1 ' / 5 , 


9 1 0 1 1 2 e 3 ) ) 


4X 


SIIFFnESS and inertia PR0PERTIES*/15X92E19o7) 
SMFKmFBS MATRlX^'V(]X9l2£lOo3n 
MASS MATr I X*/ 2E10 p 3.) ) '7 , 

( 1 VRA: ,1 C.'.i matrix^;/ ( 1 X 9 l2ElCc.3 ) ) 


1 IS 


17 UmiTY MATrUX #) 

r^-D 

lIMFTC bUiO 

MiLiROUTiNii 5TirP (A, 0) 

COf^MON/AAA/PhO^H.ABsXY 

/WnO/Z!:A,E»P,W?2»Y 
UIMCf'bloN D(6.6) 
c = ir 

'■J--ZP.2 
M 1 = 1 
M2 = 2 
.■■!3 = .3 
MA-=4 
M5 = 5 
M6 = 6 

00 lOU 1=1,6 
^ DO lOU J=l,6 
I’-^O D( I f J)«0« 

0(M1 ) =12,*F«B/A/A/A 
D !M2 tMl > *~D(M1 »M1 ) 'M-A/P'o 
D(.M1,M2)«D(H2,M1) 
D{M4iMI)s-0(M1»M1) 

D (Ml iM4) ) 
D{M2,M25«4«#E#e/A 
D(M4*M2>««0(M2»M1) 

D (M2 *M4) «0{M4»M2 ) 
D(M5tM2)«2«#e#8/A ' ' 
D(M2»M5)«D(M5»M2^ ’ 

D(M5.Ml J«0<M2» «S 

I* ( -'-V;) ’* i (': j ) 

' ■ 

D I ^ I i;, ij I p ty^ 


I- ' H » { (•'{ .1, ' ) ,k -1 ,6 ) »J = 1 ,6 ) 

ir^MAi ,( -ytiffmess matrix*/ { iox,6e:ioo7) ) 

or:.Ti''*M ■' ■ vx/x;'j'Y' ■ 


11 

111 


119 


M6 = 6 

DO ]!., 

Do 10'. J-a j6 
luO c ( [ ., 1 ) ro:.' , 

D = IJM0?'AB;:A 
C ( Ml •. Ml ) -- 1 P / 3 5 r, 
c ) r;,.: (Mi ) 

r (HA .>fil ) -0 ,)!-r/70„ 

C ( 141 .,MA. 1 (M/j ,-.'ii ) 

C «R-!<-A/210 

c ( 0 L ) --:C (m;; ,i4i ) 

C ( M‘‘} ->^4 ) =-r { I'lZ i-Ml ) 

C (M4 5f-'!&)-C:(M'4 904) 

C (M2vf-i2)=!O(-AiiA/l05c 
C { M3 ) =C (M? ,iV!2 ) 

C {,M5 ♦M2 ) =-P’»A*A/140„ 

C ( M2 sMS ) =C (M5> ♦Mp ) 

C ( M4 ♦ M2 ) =“C ( M 1 p fis ) 

C ( M 1 ♦ 5 ) ~ C ( M 3 Si (M 1 ) 

C (mSpMI) -- 13o«'P*A/420o 
C (M2pM4) =C{M4piM2 ) 
C(M3»M3)~XY 
C(M6*M6)=XY 
C (M3»M6)=XY/2<, 



C(M6»M3)=XY/2p 

ppint 11 ♦( (cn pj ) 0 = 1 ,6 )♦ 1 = 1 sS) 
li Format (40x»>fBAsic inertia ,matrix*/(ioxp6E19o7) > 

return 

END ■ i . 

SIBFTC SU03 

SUSiRODTINE’ MAT!Nf;(A»N)''||f#fli 
D 1 mens I ON' a { r 3 j S;pM| 

: Do 4 4U I • 1 ♦ N ■' '« V 
A { I » I 5 »*1 » /A ( I »J J 
Do SOU 

6ua»,9oa.,.60.ci 

fuO p- ^ 

t 'It'.'-- ..'"-i-)''' ^ ■ 

( f ■ ( . ) - 1 ) / 0 c o 02 r» 2 0 0 0 

?u VG f^ l .. I ) -- A ( < •« ,1 ) +A ( K » I ) A C I * J ) 
r.20 i/.-A,. ' 

A ( ■ Vf 


CE»RHO»YpX: 






1 20 


D I MI-NS [ ('■'N AK ( 93 » 93 ) » AM ( 93 5 93 ) 
(, f / 9 C L / AK » AM » S s » Ad 

>[ --y 

ak, ( 1 1 . 1 i ) -AK ( 11 , n )+/ 4 ,.-)^x 
AK' (I 3 H 1 ) =Ak:( n, T] )- 6 .,.-;^X /5 
A K' ( j i . T 1 ) "AK ( J 1 ? 11 ) +29 *X 

aK( J'- •, t i ) --^Ak ( J;'i . 1 1 )+6^),'X/S 

ak ( r T 2 ) =AK (I 2 , 12 )+Y 
Ak ( .J2 r. 12 ) =AK { J2 W.;: )-v 
AK ( I I •. [ :n ^AK { 1 1 , I 3 ) "6 .. +x /s 
Ak ( I : 5 » I '4 ) =--AK( 13. 13 ) + 12 .+X/S/S 

AP { J] »r-M=AK(.Jl »I5)-6o-)!-X/5 
ak { .! 3 5 I 3 ) = AK ( J3 - I3 )-12 .^fX/5/S 
Ak ( n -.Jl ) =AK ( I U Jl )+2.,-sfX 
ak: ( f3».Jl ) =AK( 13, J1 )-6o+X/S 
A ( , ) 1 « J 1 ) - A K ( .,J I t . 1 1 ) + 4 Cl '5f X 
Ak! J 3»J1 )=AK( j3,Jl)+6e+X/S 
Ak ! I2sJ2 ) =AK! !2,J2)'-Y 
AK! J2»J2)=AK( J2,J2)+Y 
Ak ( II »J3)=AK( Il»J3)+6o^fX/S 
AK! I3»J3)=AK( I 3 , J3 )‘"12 o-»X/S/S 
AK! J1»J3)=AK( Jl,J3)+6„*X/S 
AK! J3»J3 )»AK.! J3, j3)+12o*X/S/S 

cccccc 

X»ftHQ*AP^S' . , ■ ■ ■ ,■ 



Am! II# II )»AM( 11 » 1 1 }+X#S4fS/l05o ' 
AM! 13# 13. )wAM( 1 3, II I-ll e*X*S/2l0o 
AM! Jl # 1 1 J«AM! Jl# n 5-X#S##S/lA0« 
Am! j3»Il)eAM! j3»U)-13o’>fX#S/420o 
AM! J2»i:25aAMC J2#I2KXY/29 
AM! 12# IJ I«AM-!,Ilf I2I+XY 
Am! II, »I3j»AM{|l»|3)-#i.U#X*$/2lO,, 
am ! 1 3 H'3 I »»Am ! 1 3 f 1 3„l ^*13 #.*X./ 33 « 

AM ! JU I 3 ! «AM ! JU I 

AM! ] i. ».)] ^ ’•■‘••.••.•K t < ,Jt 3 -X*S*S'/ 140W^ 
I-Vf- { i 3..Ji I ) + 13,-»X}?S/420„ 

I ,.]'l ) 1 .-.JJ )+X!?S«o/lU9,, 

i'U J'l 3 + 1 1 cXX+$/2lUo 
A.’-M ! K \y 1 -.-.+ ! I . J2 3+XY/2o ; 

Am! ^j,? ,.,J2 3 =Am { ,l2 , J2 1+XY ■ 

■‘K t i !',/’} I -^ '; ! i I •-> . I '■' ) "• 1 3 1 ■>f’X 5f 3 / 42 '-3 « 

A” ( ' 4 ■. ,!3 ) "AM ( r .4 - .J3 3 +9,-.+X/7Uo 

•■■■-{ I i ' 3 A;4( I i , J3 3+1 1 „4«X+S/230o 

4 ( i “ j 4 , J 3 3 + 1 3 <, *X / 3 5 « 


i f 1 1: i-' ,i ( • L f ! ( A , |V, N 9 M 0 , E I O'f 
a{'’',!.93) ,EIG.(6) ,X!93: 
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etc A-MATRlXsMN-ORl^tRsElG-EuVo sMO-nOoOF VAi 
P I = 1 c 1 4 1 fj 9 
Do lOOn IV=1,MU 
I TMAX = 30-2^aV 


fvi=MN-I\/+l 
DO ICU I=1,M 
lUO Y(n^l. 

TTFRA--: 

l5o DO l]u J^lsN 
X ( J ) = w , 

Do ilD L-IsN 

llo X ( J) =X ( J ) +/-\( J »L ) 5!Y ( L) 

■fF( ITL'rA.GToIo) go To 12U 
M = 1 

XMAX==ABS ( X ( 1 ) ) 

DO 13D I=1,N 

I F ( XMAX o GH o AB S ( X ( 1) ) ) GO TO 130 
XhAX = ADS{ X ( I ) ) 

M=I 

130 CONriMUE 
12 0 continuf 
R»X(M) ' 

DO 14U 

140 Y{n=X(I)/R 

itera»itera+i 

IF{ ITErA^LEo ITMAX) 60 TO l5u 



EIG( IV) PR 


SS»SQRT(1«/R) 141591 

SS*100000<,#S$ 

print 17? ss ^ 

17 0 format { 4 wXf,E 19 » 7 y. r;:| 

c ^ . 0 

YMAXaY(l) ... 

M.l is .,!* 

; - DO . 

itO'-GY-XI I |*ACM?;I 5/CP^'YMA'X,) 

- . DO'' IBO, ; 'To: 

laO Ad ., nd> ( i ..))-R-;^Yd )«X( J) 
n- (Mci.fidd GO To iuOO 

filBSl'S'®# :B :> B' yV 
lldNlB-O, 1) GO TO 220 






2u 0 A ( I .>J!=A( I sj+i ) 

DO 10 I - M,M1 

DO 2 10 1 = 1,. Mi 

21 ^ Ad ,J)=A( f + l,j) 

1 FCRf.'iAT C lOX »E19o7 ) 

lu 0 COMT u;u:' . 

CiC' TO 23C 

2«^0 ] fd ' lo lQi, X ) F I X. ( i-’M ) =A ( 2 52 ) 

IF(NUl-Ot, 2 ) Elc; ^^^M)=A( I 5 I) ■ 

2% RETUPM 

emd 

bUbpooT d.'E E I ( Ak 5 Am sMm »M0 aX 5 E lG ) 

tllMiiNOlON AK{40,AG) 5AM(40,A0) 5X(4056) 5D(405A0) sY(40) 5EIG(6) 
P T = '3c 1419 
P = 0 o OOu 1 
Do 100 I 1=1, MO 
ITERA = C/ 



CA1.,L sweep CAk sAvisDaWM, I I,X) 

DO 110 i=i,mn 
llO X (I , I I ) =1« 

99 KC=0 

160 DO 120 1 = 1, mm 
Y ( T )=0. 

00 12U K«1»MN 

120 Y d )=Y( n+D( I ,K)^-XiK» I I ) 

R=Y('l) ^ : '■■Ood-: 

00 90 ■■ ■ ;.v2:,dS|a 

Y(i)«Y( n/p .. 

IF{AeS(Y( n/X(I,in-*l^O)-E) 90»9g» 
KC«1 

Continue ■ 

ItERA»i:TIiA+i:„ 

IF{ iTif?A.Ofe9 2oyi 

f P ( fC Cl 1 0 0 0 ♦ i OOfJ f t U4,^r-" 

DO d T =» 1 * MN i 

X (I 5 d ) = Y ( P, ’ 

1 \ I d '^Xat 1 ” .7) '' ■- 

--bO'.'T ( P> / !'> )/(/., »f' I ) 

isr.Dvl v,(^ (.), J d 5j=l»MM) 

2 Kr'-'-'AT dXal'O'iPdf) 

lU 0 Crq':r I .VUE 

1.1. ’r fC vSUB 11 -. 

.‘i liPPOU T i ’ i u oW E Dp ( Ak 9 Am , D , » 1 1 dX ) 


V? 


102 7 


102 

^0 


1 04 
1 1 6 

I j 00 
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D 1 mEn I ON Ak ( ) s Am ( 4^ 9 4^ ) ^ D ( 4^ » 4 U ) ? X ( 4O ? 5 ) j pQP ( 5 ? 40 ) j Sop ( £ 

I I = I I I-l ' ' " ' 

IF( II ) 111. 111.9112 ■ . 

112 DO ion 1=1 9I I 
DO 110 J = 1 .Mf.; 

ROR( I 9j)=0. 

DO lit.' K.^l?;':'! 

11 0 ROR( I 9 J) =RCR( I .J )+AH(K9J)rrX(|<,I ) 

10 0 Continue 

CALL MATINtRO^- s I I ) 

Do 130 J= 1.I I 
DO 130 I = in .MN 
SOR ( J 9 1 ) . 

DO 130 K = l» I I 

13 0 SOR( J»n=SOR( J9l )‘-ROR(Js.K)'»^ROn(K»I ) 

no 140 I^IU'^.N 
DC 14u J=1,MN 
0( I »J)s=Oo 


IFU-II) 140.140 9 141 
141 DO 150 K=l9ll 
150 D( I 9 J)=D( I »J)+A|<( I 9K)*S0p(K»J). 
D ( I , J)=D( I 9j) + AK ( I »J) 
i4o Continue 
GO TO 160 

111 DO 17U I^l.MN 


DO 170 J=1»NN 
170 D( I » J)=AK ( 1 9J) 
16 0 RETURN 
END 

c 

C 

$16 JOB 

SIBFTC TAHOEN^ 

C ' ' 

■'OD". 


c 

c 


' • ■ ■ ,rr'- t.wwstw.fiiSisi.'f**'* ’ ‘fX’' 

DImEnSi’o''' AmAS4,( .■C'{40)"9XXX(40.40) 

D I mE N iS I n*.' U ( 4'o )4 Y (,40 6 )' » E 1 ^ 6 I 

0 AI..L F LOV t O '-' ) 

CAt. L .-LDf:i-20 .'-v) ■ 11 ;:,; ' ' i 

CC'OMni.,! /!■: ANp..'i'.: /PHO»H #AB 
CON.'-iON /!■•■•■■):■• ‘.A/4:M )i E 9 R 9 022 
COM-'OH /ii!-U|. T/-.' Ar^zni 
! /A; i( >40) 

(■■ PMMfiM ./ ' ■ / \' ( 40 , 40 ) 

UmIA - 1 i f-iAX •.'■■•UO 
■MO nr. "I 


THE turBi.mE Blade 


/290i'i0'J0C o s 15 9 O 0 283 


A = 3oa4;.v 

A B 4 

Al* -■ 1 . ! c^j 12 5 

A. ! i 1 — 0 c 'J i ? 1 



] '1 1 f 


TF(w.,l'Orl) GO To 310 

cof'iTpiiir 

zp] =ZE<li'-X 

zpl = l c. 16&/12o-s«-0<, 125*^' 

W22 = W2if-W2 


Do 10'-: r-i,MN 
DO 100 J=1,MN 
A M ( I .J ) ™ 0 0 
Ak ( I » J )=0« 

S = 19,.301/25»A 
S=0c790 


H=A/AN 


Do 110 I= 1 ,M 

R*RC+H#AI 
11=2*1-1 
I 3=1 1+3 

o:all stiff (h,stif»ii 

CALL MASS (H»AMASS»II 
DO 110 J«n*I3 
DO 110 L*II,*.!3.m ,, 

AKJ J»l )»Ar,( J»li,'')'+STIF( 
AMC J*L)«AM(4»UHAHASS 


f "i 2 * ¥M^i^ '' I I i I'L 

A:■^ ( T 1 J ) AM( I + 2 > Jf.i ) ■ 
Air ( i » ,i ) "AP ( [ f 2 ) 

Do 1 3 0 , 1 = I » ■'■(14 ^ /.vj.,. ; V, 
I'G 1 Ar ,|~.l.jM:t f: 

1 '<■' ( i ! i .■ .A”^io j 



I F ( N « EO . 1 ) MH»M2 • . ■! 

f'Alt. MAT TP (AK»MM) 

ITO 7 70 ■■ ' 

P0’: 

C ( J) - C. { .1) ^Ak.( ( ,l ) »AM(L»J) 

!! ( •.■,•. pi" ) S' A V MA ) 

' '‘I..!, * 1 'V ■•;';!•■ ( »A i > V Ai'^ 9 (■';r-i 'j 6 > A X X 9 L I G ^ 
':. ,"( ! ; I I ( ; I' T* ( A r 9 aTti » 6 9 Li I G M ) 

! r;??As7 ( 2x,6ri9.7 ) 




K l2 = K2+r 


Ap;(KI2‘»XJ2)=AK( I 9 J) 

AiM( K IX 9XJ2 ) =am( I ,J) 

AL'iJlTIr.M OF GAmD'STIFFo AnD MASS MATr I CES Dqi^E UpTo SToR 
I =2’?'-M'!t'N 

DO 170 I 
1 = 1 + 2 
J = J + 2 
K =2 + M* n 
l. = K*--l 

CALL band (EsRHOjI 9 J 9 K 9 I ) 

T 1 = K 
J1 = L 
Kl=I+2 

L 1 J +2 , 

CALL BAPD(E‘ »RH09ll»Jl»Kl’Ll) 

Continue 

IF(W2oE0 6 (,)p) GO To 160 

stiffness Correction made to account Co Fo due tq mass of 

Do 320 I = 1»N '■ 0 , 

J*2#M*I“*X ^ ' ' ' '■ D-'^sfote 

AK( J»J)*AK(U»U>“29*RH0*AP#S*W2*W2 . 'Ip.l ^ 

Continue .-I 




n n rrn n 
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8 Blades = *|Iio/iox;«^bamd to blae 

11 hrrJyAT ( 30/ ;Hf rsLDUCEO MATRlCEp 

9 I'O'RMAi (E<o-X3-!fAobEMiiLED STIFF- AniL) mASS fjjATpICES-st- ) 

lU FOPMAI {/(lX.l2r.lOo3) ) ■ 

7 Forma 1 ( /ioX.^f-l-JNL)AMtNTAL MATn'^RcOo = * , 1 1 , E 19 o 7 / (1 5 X I GHEp ra' 

= '+i<F19=,7)) ' 

lA Format (/ox;.*potop speed — -sf 

3i Format ( / . luX 9 !:;Ei 3 o5 ) ) ■ ■ ! 

STOP 

END 

$IBFTC SUEil 

SUBrooTImE stiff (a»D5M) 

CO'MMOM /KAnpUr /PHOsIUAB 
C OiMMON P? OOi ! A / iB 1 9 E 9 R , IV 2 2 

dimension DCi4:>i4) 

C = E 
B = ZB1 


- length OF 'iLEMENTS 
• MOivlEN'T OF INErTiA 

mo^Ijlus 


- YOUNG 9 S 


A 
B 
C 

M - NUMBER oF ':!i EmENT 
D - stiffness M-sJRIX 
D(M»M)= 12 c,^C*B/( A*A*A) 

D(M+1».M)=D(M9(M)«-A/2o 
D ( M 9 M"*" 1 ) - D ( M + 1 9 M ) 

D ( M+ 2 ♦ M ) *-D ( M 9 M ) 

0 ( M »M +2 ) =D ( M +2 9 M 

D(M+3»M)»D(M^>I»M'5' '/.r ' ' 

D(M»M+3)*D(M+3H1) .,'.1 

D ( M+ 1 9M+i ) a44#6*C/A ■ f-Wif.'l'jY 
D ( M+3 »M+3 i «0 

0 i M+ 3 » m' I ) «D I 1 V M+ 1 /2,^i 

D iH+3 *5^+2 ri.*P'T^gf 
D _( M+ 2 » M+ 3 ^ *p { H 

■ ^ 9 M)-A^fA^f { (-P/lS. + llv / 210 . ) /A+I'^/ 105 o ) 

3 , w ) ( ( R / 10 - j 3 ^ / 42 0 o ) / A+ 1 1 o / A20 O 

M-f-i ) -X*A-;'. ( { r-V lOr+1 1 o /? 1 C ^ ) /A+3 o /70 . ) 

M t'l sM+1 )-X-i(-A-!<-A-!i •{ lo/2lO,. + i o/(105o*A) ) 

M +9 ,M,+ i ( ( - 9 / 1 U 0 +I 3 . / 420 .. ) /A- 3 ^/ 7 U j ) 
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D ( M+1 f M+2 ) =D ( M+1 ( ( 3/1 0 o + l3 „ / 420 , ) /A+3 1 ^ /420o > 

D(M+2s'M+2 )=D(M+2 »M+2)-X*( ( -R/2 « + 13 o /3 5 o ) / A-1 1 o /35 o ) 

D (M+3 »M+2 )=D{!V:+3 »M+2 { ( ■/;./ 10 «-l i « /2 10 o ) /A“23c /2l0» ) 

D ( M iiM+3 ) =0 (M 5 M+3 )-X-^tA* ( (-P t <«.13i /420 „ ) /A-2 o /35o ) 

0 (M+l 9M+3 )=D(M+1 5H+3 )-“X+A*A^t (-R/60‘'»lo/l40c: )/A-io/84o ) 

0 ( m,+ 2 »M+3 ) =0 ( V;+2 s-M+I^ )*^X^(-A.,'<. ( (B/1 0 o-i 1 0 /2 lO 0 ) /A+2 o /35 0 ) 

D ( M+3 9 M+3 ) =L)X M+'-> sM+S I'i^-Jf-A-x-A* 0 / 210 + 1 / ( 1 05 0 *A ) ) 
lO 0 return 
EmD 

SiBFTC SUB4 

subroutine mass (AjCjM) 

COMMON / KANPUR / RHOsHjAB 

DIMENSION C( I4»14) ' ■ 

C p - SPECIFIC density 

c A - length of element 

C B - AREA 

C C - MASS matrix 

B = AB 
R = RHO 

C (M>M) =l56o + R*A^(-B/420o 


C (M+29M+2)=C(M»M) 

C ( M+2 »M) =A4c-«-R*B«-A/420o 
C ( ^ M-f 2 ) = C ( M"^ 2 9 M ) 

C (M+UM) =22o^R*A^A^8/420o 
C(M?H+1)=C(M+1?K) 

C ( |vi4*3 »m*^2 ) «^C { M**^**! 

C ( 2 ^ ) «C { f N’^2 ) 

c ( M+i »M+i ) aR#e*A»*»/ s.0:s 
C ( M+ 3 » M+ 3 ) bC ( Pl+i ;»'M+1 ) , , " + 

C'< M+3 »M+l ) *f’^R^>S#A#*3'/1.404» 
C < M+ 1 » M+3 ) «C i M+3 ♦;M+.5l 
C { M+ 3 » M ) ««t3t 
C t M'+ 2 » M+ 3 )'w«C ( 

C {;M+] »MF2 ) --AC (M+3»MJ 


C (M t 
RLTUf 


M-t-3 ) -0 { ''*■'-3 »M ) 



P 


;r.IBFTC SUlSj 

subf'ot T T n;. iv 
LOiv'A’ON/H 1 /A!l 
COM:'’'. N/l'.;-' /•'•.■•; 
COMMON /O.'':!!., 
AK(<-'IC )--A'Mt-. 
Ak { I » i )=Aif ( I 
AK (Kr I )--Ar ( 1 ; 
AK. n 9 1 ' ) ■- A!-; ( I 
Ak ( L »L ) ( 1. 

AK( J» J)=AK( J 
AK(j9l. )=Ai's( J 
Ak (|. 9 .; ) =Ak ( L 
X J = r;tlO«AP + S 
AMP --K >-/''M(K 


m:: ■ ■■ - 


X i..' { 1 .. 9RH0 9 ] 93 ^KnL ^ 

(. /,u-/j.O) ■ 

[ / 9401 : _ 

“t-'. )+■•*. +2' :Oipi / S 

1 1. ) + A 0 t '•<■ Z P i / 2< 

■> I )+2:,*\-.*lpl/ S : 

» K )-t-2 o-:i-E*:!-Zpl / S 
• L)+AP*E/S_ .j- 
9j)+AP*E/S 
«L)-AP+'E/S ‘ X 

, j)_ApxE/S ■ 1 ' 

!*k)+X,j-}f'S)fS/l05“ 



Am ( I » 1) =Am ( I W )+X.H!-o#S/105 
Am( I »K ) =AM ( I SK )*Xj-)!-S*S/l40 
Am ( K » n =Am ( !< -j I ; -X ].4C.' 

am ( 1, sI. ) -A’l ( I S !. )+Xj/ 3o 
AM{ Ji.J)^--=AM( .UJ)+XJ/3„ 
AM(..j9L)=AM( Js’D+Xj/b, 

Am ( L ■' J ' =■ Aim ( L » J ) +X j/b 

RETUr^M 

FmD 




Thin book Is to be returned on the 
date last stamped. 





